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INTEGRALS AND INTEGRAL EQUATIONS IN LINEARIZED WING THEORY ‘ 

By Habvard Lomax, Max. A. Heaslet, and Fhasklts B. FunnEa 


SUMMARY 

Tki formulas of subsonic and supersonic wing theory for 
source, doublet, and vortex distributions are reviewed, and a 
systematic presentation is provided which relates these distri- 
butions to the pressure and to the vertical induced velocity in the 
plane of the wing. It is shown that care must be used in treat- 
ing the singularities involved in the analysis and that the order 
of integration is not always reversible. Concepts suggested by 
the irreversibility of order of integration are shown to be useful 
in the inversion of singular integral equations when operational 
techniques are used. A number of examples are given to 
illustrate the methods presented, attention being directed to 
supersonic flight speeds. 

INTRODUCTION 

One of the most fundamental approaches to the analjiical 
investigation of linearized ving theory, throughout the sub- 
sonic and supersonic ranges, stems from the use of certain 
elementary mathematical expressions that are identified 
uith the physical properties of sources, doublets, and ele- 
mentary horseshoe vortices. By means of these expressions 
boimdary-value problems involving vings vith thickness, 
camber, and angle of attack can be solved. These problems 
naturally fall into tvo categories: one, involving bodies 
with symmetrical thickness and no lift, is analyzed by means 
of source distributions; and the other, involving lifting plates 
without thickness, is toalyzed by means of doublet and 
vortex distributions. 

All these distributions require the treatment of singu- 
larities in the mathematical analysis. Thus, for subsonic 
Mach munbers, the concept of the generalized principal part 
pla3's an important role in the calculation of the induced 
velocities in the plane of a vortex sheet. In supersonic wing 
theory, the generalized principal part is again used in the 
analysis of vortex distributions, and it has further applica- 
tion in the treatment of conical-flow problems. However, 
the existence in supersonic flow of pressure discontinuities 
(due to Mach, or linearized shock, waves) brings about another 
type of singularity the mathematical analysis of which leads 
to the introduction of the finite-part concept. The integrals 
in both subsonic and supersonic wing theory thus require 
careful attention to the discontinuities in the integrand and, 
as an illustration, indiscriminate use of such standard devices 
as inversion of the order of integration can lead to incorrect 
results. 


When direct problems are involved, that is, when prescribed 
functions are to be integrated (as in the problem of finding 
the pressure on a wing with symmetrical thickness), a guide 
to the proper method of calculation is often furnished bj* 
physical intuition. However, when inverse problems arise, 
that is, when integral equations are to be inverted (as for 
the flat plate of arbitrary plan form), the mathematical 
methods are more abstract. Nevertheless, the solutions to 
several types of inverse aerodynamic problems have been 
obtained by reasoning that required an understanding of 
the physical nature of the flow field. This method of solu- 
tion may be sufflcient for the particular problem involved 
but it is difficult to generalize. By using the aerodynamic 
data to construct mathematical Imundary-value problems 
requiring the inversion of singular integral equations and by 
obtaining these inversions from a purely mathematical 
(operational) basis, a technique evolves whereby the existing 
solutions for two-dimensional subsonic, and three-dimen- 
sional supersonic wing problems (e. g., thin airfoil, conical 
flow, and Eward solutions) are synthesized. Furthermore, 
the solution to the general supersonic wing problem is 
suggested. 

The purpose of the present report is: First, to review the 
formulas of linearized wing theory in which source, doublet, 
and elementary horseshoe vortex distributions are introduced 
and to relate these distributions to the pressure and vertical 
induced velocity in the plane of the wing; second, to present 
an operational technique that can be used to invert the sin- 
gular integral equations appearing in the appheation of the 
above formulas; and finally, to present certain special 
examples which will illustrate the basic concepts. 

LIST OF laiPORTANT SYMBOLS 

chord of a wing 
drag coefficient / .. — \ 

pressure coefficient ^ 

complete elliptic integral of second kind, modulus i* 
complete elliptic integral of first kind, modulus k 
complete elliptic integrals of fihst and'second kinds, 
respectively, with moduli 
moduli of elliptic integrals 
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Complementary moduli (Vl— — 
lift 

Macli number in free stream 

slope of iving leading edge 

cotangent of angle between the and x axes 

cotangent of angle between the ^ and x axes 


loading coeflBczent (pressure on the lower surface 
minus pressure on the upper surface, divided bj' 
free-stream dynamic pressure) 
characteristic coordinates 

[(x-x,r+(l-Mo‘) 
wing area 

maximum thickness of a wing 
perturbation velocities in a, y, z directions, respec- 
tively 

free-stream velocity 
Cartesian coordinates 


angle of attack of wing 





A 

X 

Mi 

Ml 

M 


Vc 


’ll 2 
Po 

T 


Wc 


jump in value of the quantity considered across the 
2=0 plane 

stream wise slope of surface 

yrw ^ ... ... . 

Vl-f-m? 

1-f m^ 

1— m3 

X — X\ 

(2/-yi)*+2* 
obhque coordinates 
density of free sti;eam 
area of integration 
perturbation velocity potential 

X — Xi 

3V(y-J/i)*+2* 


SUBSCRIPTS 


I value of a quantity on the lower surface of a wdng (2=0 
plane) 

u value of a quantity on the upper surface of a iving (2=0 
plane) 


PART I—THE THREE FUNDAMENTAL FORMULAS 


normalized form, Laplace's equation in tlircc dimensions. If 
a sufficiently thin wing at a small angle of attack is situated 
on or in the immediate vicinity of the xy plane, the bound- 
ary conditions in the resulting linearized theory^ may^ be as- 
sumed specified at 2=0 and, by means of Green’s theorem 
(see, e. g., reference 1), a solution to equation (1) can bo 
written in the form 


1 ^00 p j 8*2"1 

<p{x,y, |^A«j(a:i,yi) --A»<Xi,yi) dxidyi 

( 2 ) 


where 

and 


r.= { (a:-ri)*+(i-il4»)[(y-yi)*-h2»] }'/» 


Equation (2) relates the perturbation velocity potential at a 
point (x, y, 2) in space to the discontinuities in the potential 
and vertical induced velocity at the “plane of the wing.” 
Thus, Av>=¥>ii —<pi and Aw=ta„ —Wi, where the subscripts « and 
I denote conditions on the upper and lower side of the xy 
plane. 

In a later section, equation (2) will be used to obtain ex- 
pressions for source, vortex, and doublet distributions in 
subsonic fliow. 


FIELD EQUATIOK FOR SUPERSONIC FLOW 

The form of the basic linearized partial differential equa- 
tion goveraing supersonic flow fields can be written in terms 
of the perturbation velocity potential as 


{Mo^—l)ipxt—<Pitii~'pi,— 0 . (3) 


Since M, is now greater than one, equation (3) is, in its nor- 
malized form, the wave equation. A solution to equation 
(3) that relates the potential in space to its jump Atp and the 
jump of the vertical velocity Aw across the 2=0 plane has 
also been derived by means of Green’s theorem. , A form of 
such a solution, due to Volterra (reference 2), can be written 


p(x,y,z) 


where ^ 
and 




X — Xi 

3v%-yi)*+a* 

(4) 




rc= { (x-z,)*- (A//- l)[(i/-y,)*-f 2*]}*/* 


SOME BASIC MATHEMATICAL PORMCLAS 
FIELD EQUATION FOR SUBSONIC FLOW 

The basic linearized partial differential equation governing 
a subsonic flow field is derived under, the assumption that 
perturbation velocity components are small relative to the 
free-stream velocity Vo. Written in terms of the perturba- 
tion velocity potential (p(x, y, z) the equation is 

(1— M<,*)^e*t-(-^fW+^p«=0 (1) 

where Af„ is the free-stream Mach number and the x axis is 
parallel to the free-stream direction. Equation (1) is, in its 


The area t is, that part of the 2=0 piano contained within 
the Mach forocone from the point (x, y, 2); that is, the 
area bounded by the line Xi= — <» and the hyperbola 
(x-zi)*-3*(j/— yO*— j9’z*= 0. 

Obviously equations (2) and (4) are not similar although 
the basic equations to which they apply are. A foimal 
similarity can be obtained, how'over, through the intro- 
duction of an integral operator, originated by Hadamard 
(reference 3), and referred to as the “finite part.” The 
use of the finite part and certain other techniques necessary 
to reduce equation (4) to a foim similar to equation (2) 
requires some attention and a discussion of the mathematical 
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difficulties iuTolved vdll be giTen in. tbe foUovring section. 
Tben an appKcation of these techniques in subsequent 
sections will make it possible to obtain expressions for 
‘source, vortex, and, doublet distributions in supersonic flow. 

THE FINITE PART OF AN INTEGRAL 

In the study of linearized supersonic flow problems, one 
is continually confronted with expressions of the form 



The integrals are of this form in the sense that the integrand 
is infinite at one (or both) of the limits and this limit is a 
function of the variable by which the partial derivative 
is to be taken. Such an expression is annoying because 
the derivative cannot be “moved” through the integral 
sign according to the usual rule, namely, 





and it follows that 


r A{y)dy 

X (x-yr^ 


X 


‘^ A(x)—A{y) 

(x-yr^ 


dy+A(x) 



dy 

{x-yf‘^ 


( 10 ) 


The final generalization of the definitions given by equa- 
tions (7) and (9) is accomplished by considering the 7^ de- 
rivative of the integrals and, furthermore, by allowing a 
functional dependence on x of the integrand A. First con- 
sider the definition 


a V ACy)dy 
jXvSa:/ -yjx—y 


<-!)■ 

\^xj X -yfx^ 



A( 3 /)dy 

( 2 ;_y)n+l /2 

( 11 ) 




F(x,y)dy=F(x,x)+ 


^F(x,y) 

bx 


dy (6) 


Direct application of equation (6) to equation (5) obviously 
yields an unacceptable indeterminate form since the term 
corresponding to F(^x, x) is infinite. One way of avoiding 
the difficulty is to integrate equation (5) by parts so that 
the radical appears in the numerator of the integral and 
then to apply equation (6) to tiie resulting expression. 
Such a procedure can be carried out without the iutro- 
duction of any new mathematical symbol or concept. 
However, this involves unnecessary restrictions on the 
integrand and often leads to unwieldy forms since the 
derivative of the funtion f(x, y) with respect to y can be 
cumbersome. 

Definition. — more direct way of taking the derivative 
throu^ the integral sign in equation (5) is accomplished by 
using the integral operator known as the finite part. Con- 
sider the simple equality 


A p _ 1 

-y/x — y -yjx — a 

The finite-part sign can be introduced by the definition 

dy C‘ dy 

X y/x—y SarJ„ -^x—y 


(7) 


from which it follows 


f 


dy 


—2 


(®— y)®^* ijx—a 

The natural extension of this idea is to consider 

b A( 3 /)dy _ b A(y)dy 
X ^x yjx—y J<i -\ix-y 


( 8 ) 


(9) 


where A(y) is continuous at x=y and is integrable elsewhere 
in the range of integration. The evaluation of J can be re- 
duced to a form that requires only the definition introduced 
by equation (7). Thus, by adding and subtracting the 
same term, J can be written 


The second integral can be expressed in the form 

p Ajy)dy r‘’ A( 3 /)—B(x,y) j B(x,y)dy 

X (x-yr^^^^~X (x-y)’+^^ “^"^X Cx-y)"+^'^ 


where 


( 12 ) 


and 


-B(r,y)=A(x)— A'(ar)(a;— 2/)-(- . . . + 

r dy (- 1 )* 2 * /'byr dy 
X (*-y)*+^'*~l-3----(2i-l) Wi 
-2 

(2i-l)(x-a)*-'/* 

Finally replacing A( 2 /) hj A{x, y), equation (11) again defines 
uniquely a finite-part integral provided that 


Methods of evaluation. — ^If expressions of the t 3 rpe pre- 
sented in equation (5) appear in an analytical development, 
it is now possible, by using the finite-part symbol, to take 
the partial. derivative operation through the integral sign. 
Such a process needs no further amplification. In applying 
the results of such an analysis to the solution of some spe- 
cific problem, however, one is confronted with the inverse 
operation, that is, the problem of evaluating the finite-part 
integrals. This can always be done, of course, by means of 
the definitions already given. Often, tho ugh , such evalua- 
tions can be simplified by using one of the two following 
processes. 

The first process is readily outlined. Eewrite equations 
(8) and (10) in the form 

r^ J(x,y)dy y f(.x^y) — /(z, x) , 2j(x, x) 

X (^-y)®'®“^j. (X-yr^ 

and set the indefinite integral of 

M{x,y)dy 
J {x-yW 
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equal to F (x, y) + C. It foUoivs that 

I 

■[F{x,a)+C] (13) 




whei’e 


(7=lim \^^-F{x, j/)1 (14) 

v^L-sx—y J 


Tlie secood technique avoids the necessity of evaluating 
the constant C. It depends on the use of the complex vari- 
able and is valid only when /(x, y) is real in the interval 
a<y<6 where h is some number greater than x. Again set 
the indefinite integral of 




jx,y)dy 
(x—i 


equal to F (x, y) + C. Now if r. p. stands for the real part 
of a function, the evaluation of the finite-part integral is 
provided by the equality 

( 15 ) 

As an example of the second technique, consider the 
integral 


I=-P 


y^dy 


. „ _£‘ y^dy 

■ J. (x'-j/*)’/* {x^-yr<^ 

where a:<5. From the relation 
y^dy y 


fi 


• V 

arc sm — 
a: 


together with equation (15), it follows that 

b 


/ b . b\ v 

A simple extension of this result yields 

r- fdy p T* 


.y2)8/S 


£ 


(fio+Ciy)dy 

yf7i 


■j:. 


(Cp+Ciy)dy_ 




-2 


q -yja+bt+ci^ 
Included in this result is the equality 


£ 




jrs/i 


dy=Q 


which contains the very important identity 

dy dy ^ 

Xx, (x'a-j/i^/Vy-J^i Xi (l/— 1/ 


The case of multiple integrals. — WTien applied to the 
analysis of single integrals, the above defintions of the finilo 
part cgincide with, or are a ro-expression of, those given 
originall}" by Hadamard (reference 3). Hence, 


rfix, y)dy r.f(x, y)dy 

J. {x-yy*~Ja i.x-yf>* 


ivhere | was the s 3 unbol used bj' Hadamard to denote eval- 
uation b}’- finite-part methods. When applied to double 

integrals, however, the signs f” and are no longer equiv- 
alent. Hadamard, as well as A. Robinson (reference 4), 
maintains the convention that the order of integration in the 
operation | is reversible; that is,* 


IJ dyjdxj{x,y)=\^^dx ^dy J{x,y) 

Such a convention requires that all singularities for which the 
j order of integration is irreversible must be excluded from the 
area of integration. These singular regions are then treated 
separate!}'. This convention has the disadvantage that, in 
evaluating multiple integrals, the value of a given integral is 

not independent of succeeding integrals. The operator -j- 

avoids difficulties of this kind and each definite inlegnil is 
independent of succeeding operations. For example, 

n dn _ dy _Q 

X J' Jo (?— 

but, according to reference 3, 


An application. — The above methods can be applied to give 
the following simple but useful result. Let F =a-l-by+cy‘= 
(Xj— y)(j/— X])(— c) and q=4a^—b^, then 


[Ctf(2ct-j-6)— Ci(6t-H2a)] (16) 


(17) 


( 18 ) 


f — 

Jo a- 


dr) 


= 0 


(?— 

whereas, accoi'ding to the same reference, 


de 

Jq Jo 


dv 


(l-^)’^S^ 


‘2t 


Although the use of the symbol makes each integration 
independent of subsequent integrations, the order of into- 

cannol be re- 


f 


gration for operations involving the sign 
versed. Hence, 

dxj{x,y)^ J dx f dy ji.x, y) 

For example, the relation 


=0 


holds while the same integral taken over the same area but 
in reversed order is 


•since the order of Integration plays an Important role In the fonowitig dCTClopmcnt, inte- 
gration flrsfwith respect to i and then with respect to y will bo denoted fdyf dif (r, y) while 
Integration first with respect to y and then with respect tox will be denoted y dry dy/Cr, y). 
When the notation f y) dydr Is used, the order of Integration b Immaterial. 
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Cdv-P 

Jo J, Cf— 




^=-2-1 


con- 


’i 

The operation defined by the symb 9 l ^ will be used 

sistently throughout the present report. Hence' attention 
must always be paid to the order but not to the multiplicity 
of integrations. It can be seen that the finite part of a con- 
ventional-type integral coincides Avith the value found by 
standard methods. 

THE GEKERAIJZEO PRINCIPAL PAST OF AN INTEGBAL 

Another type of important operation appearing in the 
development of both subsonic and supersonic wing theory 
appears implicity in the expression 

where /(y, z) and its derivatives are bounded and continuous 
in the interval a^y^b. In the attempt to simplify by 
letting z approach zero before performing the integration, a 
second special integral operator will be introduced. 

To simplify Jo, first integrate by parts 

Io= lim ^ arc tan - — -—j(h,z) arc tan - — ^ + 

*-to+ L ^ 2 

Then since 

lim f arc tan dy^=vf,{y,o)— 


U [/r(ff,o)+/r(6,o)l 


In becomes 


y—b y—a J. y—yi ‘ '' 

Definition of the generalized principal part. — The expres- 
sion of lo given in equation (20) contains the integral of the 
function/'fyi) / (y-yO . Such an integration is not, in general, 
convergent; however, when the integral is so iviitten without 
further qualification it is generally accepted that the singu- 
larity occurring in the integrand is to be treated using 
Cauchfya principal part. Evaluation by such a method is 

often indicated by the symbol ^ and is defined by the equa- 
tion 

X Vi—y U« Vi—y J,+. yi— y J ^ ^ 

or, altemativdy. 

To assure the convergence of the right-hand side of equations 
(21) and (22) it is sufficient but not necessary to assume 
that A(y,) is differentiable at the point yi=y and that 


elsewhere within the region of integration A(j/i) is either 
continuous or possesses integrable singularities. The con- 
cept of the Cauchy principal part is so well known that the 
symbol on tbia integral is often omitted, as shall be done 
here. 

The differential operator in equation (22) lends itself 
readily to a generalization of the principal-part concept. 
Thus, for the next higher order, the definition (see also, 
in this coimection, reference 5) 

(23) 

applies and, in general. 

The generalized principal part can also be expressed in 
rile form 

Jl* A(yO£fyi _ A(yO-Hfy,yO , Sj B{y,y^dy^ 

T. (yi-y)”^^~J. (yy-y^^^ (y.-y)"^' 

where 

H(y,yt) = A(y)+^(yi-y)+ . . . (Vi-y^-^ 


and 


_T» _^yi__l r Jyj__ , 

X (yi— y)*'*’'^ iiVby/ Ja yi—y 

=ir_j 

i L(5— y)‘ (a 


-y)‘} 




The first n derivatives of A(yi) are assumed to exist and be 
single valued at yi=y while elsewhere in the range of inte- 
gration A(yi) may possess integrable singularities. This 
definition is in a form that involves no extension beyond 
the concept of Cauchy’s principal part. 

It is possible to extend the definition contained in equation 
(24) to include a functional dependency on y in the numerator 
of the integrands. Thus, replacing A(yi) by A(y, yO, 
equation (24) again defines uniquely a principal-part inte- 
gral provided the first n derivatives of A(y, yi) with respect 
to y and yi exist at yi=y. 

Method of evaluation. — Operations involving the symbol 

c 

^ can always be performed by means of the definitions 

just given. However, another method can be used, which 
is often simpler to apply. If the indefinite integral of 
A( 3 fl)[{,yi—y)’''^^ exists such that 


J 


A(yi) tfyi 

(j/i— 2f) 


a+I 


■-0(3h,y)+0 


(25) 


then the value of the generalized principal part can be found 
by following the conventional rules for substitution of 
limits; thus 
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-A(yi) dyi 

(2/1— y)"'*'' 


= 0 {b,y)-G{a,y), yy^a, 


b 


(26) 


while the relation between the differential areas, as deter- 
mined from the Jacobian of the transformations, is 


The proof of this result can be obtained by mathematical 
induction. 

An application. — RetiuTung to equation (20), since an in- 
tegration by parts yields tlie relation 

X" o) dyi _ f{b, o) /(g, o) p /i>i(2/i.o)<^2/i 
X (^i-yy y—b y-a J. y—yi 

the limiting process symbolized by Ig can be expressed as 
an integral that contains only the function evaluated at 
2=0 and not its derivative. Thus, finally 




Pi/t* 


(30) 


: The value of rc, as defined under equation (4), becomes 




— , 2(i?-7;i)(f— ^i)(mima-hjS») 


/*! 


yifi 2 




/*2 


-;9*z 


- 11/2 

■] 


(31a) 


where /3®=Mo*— 1. If the variable p, is introduceil such that 


lim 

1^0+ 


5 r zjjyx, 2) 

L ^zLzi-\-{y—y^) 



/(2/i, o)dyi 
J. (2/1-2/)* 


THE OBUQUE COORDINATE SYSTEM 


^My, o) 

(27) 


Equations (1) through (4) gave the basic partial differen- 
tial equations of wing theory, together with their solutions, 
in terms of the usual Cartesian coordinate system, the x 
axis extending in the direction of the undisturbed flow and 
the y and 2 axes oriented normal to this direction in such a 
way that boundary conditions for the wing can be specified 
in the 2=0 plane. In tlie study of supersonic flow fields 
it is at times mathematically convenient to introduce in 
the 2=0 plane new coordinate axes making arbitrary angles 
with the X and y axes. 

The general case. Consider the |, ij, 2 coordinate system 
(fig. 1) such that the 2 axis is normal to the plane supporting 



Fioure I.— Coordinate systems. 

the botmdary conditions wbile | and i; are noimal to z and 
coincident with the lines x= —mgy and x=miy. If 


Mi=Vl+mi*, (28) 

the equations relating the two systems of coordinates are 




All 


Al2 


f= 


j 

fit yt 


2 = 2 , 




fit(x—miy) 

mi+mj 

Hiix + rriiy) 
wii-f-mj 


2 = 2 


(29) 


X — x^ 

"“^(2/-yi)*+z* 

it follows that in the transformed coordinates this 
becomes 


(31b) 


variable 


(i?-:;0^+(f-fi)^ 

All /l2 


rrit 


Pc = 


[(u-^0~(i-fi)£]+2* 

Finally, the differential operator b/bi transforms to 


(31c) 


5 1 / A, 

bx mi+ ms V * br\) 


(3 Id) 


The area r over which the integration of equation (4) is to 
be taken is stUl, of course, bounded by the hyperbola 
and the line »i=— oo. The asymptotes to the hyperbola 
become', however. 


17 — i?i = ^^ (£— $i) 
V-i 

)J — — (£— fi) 

Ms 


mr±i' 

m,-i8 
ma— ;8 
771 i+j 8. 


(32) 


Figure 2. shows how the area t in the xy plane transforms to 
the $17 plane. 

(Although the $ and 17 axes are oblique nith respect to the 
original axial s^'stem, no inconsistency results if the equations 
in the transformed variables arc plotted relative to orthogonal 
axes.) In case both mi and ma are less than /3 (the case for 
which the sketch was dra^ra), the asjmptolea are straight 
lines with positive slopes, and the area is bounded in both the 
$ and 1? directions by the maximum points $a, 17* and tit, 
respectively. These maximum points arc shoirai in the figure 
and their analytical e.xpressions arc 


. j. 

* ' mi+ma ’ 

t . 2Aia(mima-bjS^ 

” (mi-kmalV/S^— ma* 


Vn=n 


Vii=V 


2/xi(mimaH-j3*) 
(mi-t-ma)V;8*— mi* 

2/iiV/ 3‘— ma* 


K33) 


mi+mt 


■ If the oblique coordinate s^’stem is chosen such that mi is 
less than and ma is greater than jS, one of the asymptotes in 
equation (32) has a slope of opposite sign and the area t is 
unbounded in the 171 direction although $*,17, is still the point 
which deteimines the farthest extent of t in the $1 direction. 




Ca) xi, i/i plane. 

(b) &, in plane. 

Fhiube 2.— Area of integration, r. 


On the other hand, if m2 is legs than and mi is greater than 
/3, T extends infinitely in the ^ direction and represents 
its upper bound in terms of %. Finally, if both mi and m^ 
are greater tban the area t is not bounded for either 
negative or positive £ or 171. 

The characteristic coordinate system. — ^Another special 
case of the f,i7,2S3’stem (thes, y, z system is, of course, a special 
case obtained when mi=0 and m2=“) that is important 
enough to receive a particular notation is the one obtained 
when the ? and ij axes are parallel to the Mach lines (i. e., the 
traces of the characteristic cones) in the 2=0 plane. These 
axes are shown in figure 3, will be designated r and s, and are 
given by equations (29) whenmi=m2=|3 and fti=fi2=fi^o; thus 
213837 — 53-^ — 81 


W,_o=^o=;gV(r-J'i)(s-Si) (35) 


and the area r, shown in figure 4, is bounded by the straight 
lines Ti=r, Si=s and ri=Si= — . 



n . 

FmuRE 4. — ^Tntesratioa area Eo. characteriatio srstem. for 
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THE THREE FUNDAMENTAL FORMULAS IN SUBSONIC FLOW 

The parallelism between the basic foimulas in subsonic and 
supersonic \ving theory is so obvious that it is advantageous 
to present firet the somewhat more classical results applying 
to the purely subsonic regime. The immediate objective is 
therefore to present as briefly as possible the expressions for 
the perturbation velocity potential due to a distribution of 
sources, doublets, or elementary horseshoe vortices and then, 
by means of these expressions, to relate the pressure and 
vertical induced velocity in the 2=0 plane to the weight of 
these distributions. 

THE PERTURBATION POTENTIAL AT A POINT IN SPACE 

Tire linearized form of the perturbation velocity potential 
due to a unit source, elementary horseshoe vortex, or doir- 
blet situated in a free stream moving at a uniform subsonic 
velocity Vo is given as follows: 

Unit source ■ ; ^= — l/4a-re 

Unit elementary horseshoe vortex. tp=^zvJ4iTro 
Unit doublet 

where Ve is defined in equation (31b) and Ma- 

li is well known that a distribution of sources in the s=0 
plane splits the streamlines and forms a field sj^mmetrical in 
u, V, and <p above and below the source plane. Hence, the 
strength of the sources is related to the term Are (which, in 
turn, is related to the gradient of thickness of the simulated 
body) while the variables w, v, and are continuous. On the 
other hand, a distribution of elementary vortices or doublets 
causes a discontinuity in the stroamwise induced velocity 
(or, what amormts to the same thing, the perturbation 
potential) across the reference plane but, at the same time, 
causes no division of the sti'eamlines. The strengths of the 
vortices and doublets are therefore related to the terms A« and 
A^3, respectively, (which, in turn, are related to the wing 
loading), and produce no discontinuities in w. The exact 
analytical form of thfese distributions can be obtained readily 
from equation (2). 

The source distribution. — ^The velocity potential induced 
1>3' a distribution of sources over the 2=0 plane follows im- 
mediately from equation (2) since, by symmetry, Ap must 
be zero. In practice, the area over which the sources are 
distributed is limited to the area S defined by the plan form. 
Hence, 

<p{x,y,z)=^ (36) 


The elementary- vortex distribution. — Equation (1) was 
written in terms of the perturbation potential (p. It could, 
however, after differentiation have been expressed in terms 
of any one of the induced velocity components and the so- 
lution in equation (2) would then also be expressed in terms 
of the particular velocity component chosen. Consider such 
a case, taking for the dependent variable the streamwise 
perturbation velocity instead of Equation (2) then be- 
comes ...... 


“<*■»' *>-47/-., r.( 


bAv, 
bz To 



If the field is to be without sources, both Av) and bAm/bx 
vanish. But for an irrotational field the equality' bAw/bx— 
bAujbz holds and the first term in tlui integrand of the above 
equation is zero. B.y definition 


J_ u(xuy,z)dxi=ip(z,y,z) 

from which it follows that if the operator j dx is applied 

*J -» 

to both sides of the resulting equation, (he relation 




follows where the area of integration is limited to tlic wing 
plan form. This result expresses the perturbation velocity 
potential due to a distribution of elementary horseshoe 
vortices over a wing plan form in the z=0 plane. 

The doublet distribution. — The solution for a doublet dis- 
tribution, just as in the case of the sources, follows ini- 
mediatelj’’ from equation (2). Since the streamlines arc not 
divided by the doublets the term containing Aw vanishes. 
The doublet distribution exists, however, not only over the 
whig area but also over the vortex wake streaming down- 
stream behind the wing since the discontinuity in the po- 
tential peraists in this region. Designating the wake area 
by Tl’', the final expression for the perturbation potential 
associated with the doublet sheet becomes 



REDUCTION TO THE PLANE OP THE WING 

The aerodynamicist is usually interested in the forces on 
the surface of the wing itself and, as a consequence, it is 
pertinent to consider each of the above foimulas in the 
limiting .case as z approaches zero. An explicit expression 
of these results is given below. 

The source distribution. — The limiting value of equation 
(36) as 2 approaches zero is obtained immediately' by simply 
setting 2 equal to zero. The resulting expression is 


<p(.x, y, 0 )=^ U dxi dyi (39) 

8 

where 

?’o= [(k — a:i)*+ |3^(2/- y 1 )^'" 

Practical interest is usually concentrated on the relation 
betw'een the pressure on the wing surface and the wing shape. 
Since in linearized theory' pressure coefficient Cp and wing 
slope Xb of the upper surface are known to be 


Cp=- 


-2u _v)u 

TT' Vo' 


1 Aw 
“2 Fo 


- (40) 


it follows, after differentiation of equation (39) W'ith respect 
to X, that pressure coefficient is 



(41) 
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Tie elementary-Tortex distribulioii. — When the sti*ength 
of the elementary horseshoe vortices is known over a vring 
plan fomi, there is no difficulty in finding the potential in the 
r =0 plane since it follows from a direct integration of the 
vortex strength. The pertinent question is, rather, to 
determine the vertical induced velocity in the plane of the 
wing from the given vortex strength. If load coefficient is 
defined in the usual way 


Ap 2Au 4Uu 


(42) 


the answer to this question requires the evaluation of the 
following limi ting process: 



z(x—xQ 


[(2/— 2/i)*+z^''c g 


^ dxs_dyi 


If, as in equation (27), the generalized principal part is 
introduced, the required expression becomes 


w _ — 1 XX. (J— a?i) 

Fa 8x TT (y—yiyro 


^ dx, dyi 


(43) 


The doublet distribution. — In the case of doublet distribu- 
tions, the relevant problem requires the expression of vertical 
induced velocity in the plane of the wing as a function of the 
doublet strength A^. If equation (38) is differentiated 
with respect tp z and z is then set eqiml to zero, one finds 
without difficulty the final formula 


w 



(44) 


THE THREE FUNDAftlENTAL FORMULAS IN SUPERSONIC 

FLOW 

The purpose of this section is to repeat for supei-sonic wing 
theory the developments presented in the preceding section 
for subsonic theory. In order to maintain the formal analogy, 
it is necessary to introduce the concept of the finite part. 
This latter concept, in turn, introduces into the analysis 
integral expressions containing certain inherent eingulariiies. 
Such singularities are, by definition, points across which the 
order of integration cannot be reversed.® The study of these 
singularities and their effect on the fundamental formulas 
suggests the introduction of an oblique coordinate system 
defined in the first section of this report. Hence, the follow- 
ing analysis wiU be presented in the f, ij, s S 3 'stem, while 
transformation to the Cartesian and characteristic systems, 
it will be remembered, can be made bj' considering the special 
cases 

For Let \ 

X, y, z mi=0, « V (45) 

r, s, z mi=mi=p ) 

Following the development of the basic formulas, a summary 
of results wiU be given in terms of the x, y and r, a coordinates. 


THE PERTURBATION POTENTIAL AT A POINT IN SPACE 


The velocity potential due to a unit source, a unit elemen- 
tary horseshoe vortex, and a unit doublet is given as follows: 


Unit source <p=—l/2trrc 

Unit elementarj' horseshoe 

vortex <P=— ^Vcf2 tiTc 

Unit doublet ip=z0^l2irr^ 






(46) 


where and are defined in equations (31) and ^=Mo~ — 1. 
In this notation, the only difference between these expres- 
sions nnd the corresponding ones for subsonic flow is the 
factor 2. A. much more important difference exists, how- 
ever, in the eft’ect of the change of sign of 1 — Mq® on and Ve- 

The first task is the expression of the perturbation ve- 
locity potential at a point in space m terms of distributions 
of the elementarj’ solutions and this can be accomplished by 
an appropriate anah'sis of equation (4). Just as in the 
subsonic case, the som’ce strength is given by Aw, the vortex 
strength bj' Au, and the doublet sti-ength b\- Aip. 

The source" distribution. — The potential induced by a 
source distribution canbe obtained from equation(4) by setting 
the term containing Atp equal to zero (i. e., bj* removing all 
the circulation from the flow field). Bi' means of the nota- 
tion 


x—x. 


a>c= 


^(n-^.)+^^(f-fi) 

Ml Mg 


J3 -V (3/— l/i)*+2* 




m 


The remaining expression can be- written in the i?, z co- 
ordinate s^-stem simply as 

^(f, V, g) = 2 ~ I JJ i?i) arc cosh oiedficf 171 

(48) 


It is possible to “move” the partial derivative operations 
through the integral signs if the equalitj- 


5 

JiCt) 


/(a-, y) dy=/(a,i/)^— /( 6 ,y)||+J^ 


Mix,y) 

bx 


dy 

(49) 


K used. Obviously, if the operation indicated by equation 
(49) is to be applied to equation (48), the order of integration 
and the limits on the latter equation must be specified. 

Consider the case when the | and 17 axes are chosen so 
that both JUi and m* are less than jS. The area t is then the 
one shown in figure 2 . Further consider the case when the 
first integration is made ^vith respect to 171 . Then, equation 
(48) becomes 

J**"dfiJ^ di 7 iAw(^, 17 O arc cosh oj* (50) 

where Xi and Xj are roots of the quadratic expression r/ ( 171 ) =0 
and ^a, Va are defined in equation (33). It can be shown 
that when £i= 5 a, the equahtj^ \i=X 2 applies, and, further, 
that the integrand of equation (50) does not vanish. Hence, 


> More precisely, if an tniiereat sfnsalsrity exists in t&e area over which a doable integral ts 
to be evaloated, the difference between the values of the integration made In one onler and 
then over the same area but with reversed order Is not zero. 
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the two partial derivatives can be moved through the first 
integral sign without the appearance of the additional tenn 
involving the derivative of the upper limit. It is also not 
diflBlcult to show the value of We ac nh equal to either Xi or 
Xj is just 1. Hence, the two partial derivatives can also be 
moved through the second integral sign without the appear- 
ance of additional tenns. Finally, since 


Once more consider the case when mi and mj are both 
less than /3, the area r is the same as that shown in figure 
2, and the first integration is made witli respect to rji. Then 
equation (53) becomes 








the velocity potential at a point ($, rj, z) due to a distribution 
of sources in the 2=0 plane can be written 




mi+m. rC Awiiu rg) 

2ita(iMsJJ rc 


d^idrii 


(51) 


The elementary-vortex distribution. — The potential in- 
duced by a distribution of elementary horseshoe vortices 
can be derived in a manner analogous to the derivation for 
the subsonic case. Thus, the solution given by equation 
(4) is written for the induced velocity u rather than the 
velocity'- potential.^ Since the flow field contains no sources, 
bAufdz is zero (bj'- the same aigument presented for the 
subsonic derivation) and the solution .can be expressed in 
the form 

^ JJ 17 dxydy, 

T 

However, by definition 



so the relation becomes 


Finally, in terms of the rj, z coordinate system the equa- 
tion for the velocity potential due . to a disti'ibution of 
elementary hoi’seshoe vortices in the 2=0 plane can be 
written 




JJ^;Ati(fi,,7i)d?id,i (52) 


The doublet distribution. — The potential induced by a 
sheet of doublets can be obtained from equation (4) by 
setting the term involving Aw equal to zero (i. e., by re- 
moving all sources fi-om the flow field). Expressing the 
result in the f, rj, z coordinates, one has 


T 

Again it is necessary to carry the partial derivative opera- 
tions through the two integral signs. And again, just as 
was the case for the somces, this requires that the order 
of integration and the limits on the integrals be specified. 

* The feet that equation (4) can be written for u as well as «> does not follow Immediately 
In the case of sapersonlc flow because of the presence of the discontinuities In « ^ong Mach 
wares emanating from supersonic edges. Proof of the vaJldlty for cases of Interest herein 
is given In reference 6. 


More caution is necessary in moving the derivatives through 
the integrals than was required in the study of the source 
case, since at tlie integral is indeterminate. It is 
true that the interval of integration Xs, Xi is zero, but it is 
also true that the integrand is infinite. The value of such 
an .indeterminate form must be obtained by some process 
such as the following. 

The upper limit to the integral, fa, is a function of the two 
variables. 2 and f (see equation (33)); tha t is, in function^ 
notation fa=f«(f, s). Replace 2 in by V2*+e*, then when 
e is zero f* has not changed. But if fa(S, z) is replaced by f.(f, 
Vz*+s“) in equation (54) and the limit taken as e ajyproaehes 
zero, the indeterminate form mentioned above can be eval- 
uated. Hence, consider 


P (f , 17 , 2 ) = lim ^ 

\jU? Ot? of/ 

T , j;* 

<^li . A<!>(fi, 1 J 1 ) (55) 

By ap^ying the operation described by equation (49) and 
l etting a primed function symbolize its value at f = $»({, 
V2*+e*), the expression 






s /-?*' { 7 ,^+ 7 , ri)/.? f; «" 


is obtained. 


(56) 


The first term in equation (56) can bo greatly simplified by 
use of the mean value theorem. In the limit as e goes to zei-o 
both X/ and Xj' approach the common value of rja- Thus, 
for € very small the variation of A^' and v/ in the range 
X/<(i;i<(X 2 ' is slight. The same cannot, of course, be said of 
1/re', since X/ and X»' are the roots of re'=0. Using the 
functional notation r/=v/{rh) and applying the mean value 
theorem, one can \vrite for the fii-st term in equation (56) 


‘is ^7, 

Now from the definition of rc given by equation (31a) 

Hi p*' drij 

Jxi' fe' Jj'i' V(^-"^)"(ui-Xi0 

which is independent of e. Hence, tJio expression given in 
equation (57) reduces to 


2 A^(£i, ijo) 
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Moving tihe partial derivatives through the second integral 
sign in the last term of equation (56) can be accomplished 
by introducing the finite-part sign defined previously. 
Since 

n A+i. 

mi-f-m* \jj.2 Sf/ Tc- Te® 




it follows that 

a . 1 a 


■£ 


(gi, 

r/ \ fiiiii ) 

By means of these equalities, equation (54) can finally be 
written in the form 


For 0<mi<ft |3<mj<oo 


1a/,. ^ zi3*(m: 


(60c) 

J Jr , 


J Jr r, 

For0<mi< oo, 0>mn<i8 . 


Ay(^,^?0 


Ag>(gi,i?i) 


1 A /► \ z^*(mi+m*) r J S jy 

*’-2 ~ 


zj8*(mi+m^ 




2iTMiMa 

For j8<mi< “,/3<m2<«= 


r J ApdijTjo 

— n 


(60d) 


(60e) 


(60f) 


^ (fj n, ^)=2 


zp%7rii+m^ 
2ir[iifJL2 




Aytfi, ^?i) 


(58) 


K the area of integration is not changed but if the order of 
integration is reversed, it can be shown by a process identical 
to the one just described that 


^C^,v, 2)=o A<o'(&j Vi)- 


2irfiifi2 




A¥>(fi, ’ll) 


(59) 


The area r used in equations (58) and (59) has been defined 
as that shown in figure 2, the axes | and tj both chosen so as 
to lie outside the ilach cone from the origin of the Cartesian 
system. As the f and jj axes approach the Mach lines in the 
X, y plane, that is, as m* and mi approach P, the residual terms 
in equations (58) and (59) approach (K) A«j(J:o, — oo) and 
(X)A¥>( — <=, ij»), respectively, which represent the jmnp in 
potential infinitely far ahead of and to one side of the point 
(f, 1), z) at which the potential is being measured. Hence, in 
aerodynamic applications, the (K)^v’(fo, — °=>) and (KlA^j 
( — 03, ijj) can be taken as zero. Thus, when the 17 axes 
He along the Mach lines, thereby becoming the r, s axes of 
equations (34), the expressions for tp are without the res- 
idue terms and the order of integration is immaterial. 
When mi and mi are greater than j3 the same is true e., the 
terms 00Ap(f« v<d and (H) A¥)($s, vt) are mi s sing from equations 
(58) and (59) ,respectively) so that the effect of a distribution 
of doublets on the velocity potential can be summarized as 
follows: 


For 0 <mi<j3, 0 <mi<i3 


v=\ A(p(£„ 

(l»> 17 s) 


zj3-(ml-^-ma)XJ^ C j A^o(fi,iji) 

2irni[i2 J Jr rc 

J ‘X- J’y 


(60a) 

(60b) 


2xAtiM2 JJ V 


(60g) 


There exists an interesting corollary obtained by sub- 
tracting equation (60a) from (60b) ; namely, that the differ- 
ence between an integration of supersonic doublets made 
first in one order and then in the reverse order is equal to the 
difference in the magnitude of the distribution at two points 
in the^ plane. 

HEDCrCnON TO THE PLANE OF THE WING 

The next problem is to consider the above formulas in the 
limiting case as z approaches zero. 

The source disiiibutaon. — ^The potential in the plane of the 
disturbing source sheets follows immediately from equation 
(51) by simply setting z equal to zero; in this way 


.. ff d,, (61) 


2TCfli[l2 JJ '’0 


In order to relate the pressure coefficient Cj, to the slope 
\u of the upper surface of the wing (where both Op and Xu are 
defined in equation (40)) the operator 

2 d 2fii[i2 ^ ^ ^ 1 ^ ^ ^ 

~Va dx~ Vo (mid-TTiz) \iti Pi 2>ij/ 


must be applied to both sides of equation (61). Hence, 




” (Cli Vl) 
J’o 


dliidrii 


(62) 


The task of moving the partial derivatives through the 
two integral signs presents a problem identical to the one 
studied in reducing equation (53) to equations (60). Two 
inherent singularities in the area occur at the points £,, ita. 
and fj, 176 and form, just as in equations (60), certain resid- 
uals there. In the present discussion, interest is confined to 
the case when z=0. Equation (33) shows immediately 
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that the values of rja and fj>, rji for z equal to zero' can bo 
written 

C’?b)»-o=(’7»)i-o=’? 

It can be shown, tlierefore (the details being omitted since 
they are precisely the same as those jdescribed in the reduc- 
tion of equation (53)), that the following relations hold; 



For 0<mi</3, 0<mi</3 

3 (^J V) ~ 

V;8 — m,* 


2 


1 Ts — (ft) ’ll) 


ViS*— mj* 


v)- 


2 (ffliH-ffla) 
VH1H2 




?*0 


^(> 1 -^ 1 )+^ (?-S.) 

£i y-i ^ 

r«*> 


(64a) 


(?i, i?i) 

(64b) 


For fl<mi<j3, /3<m2< 00 


0,= 






VUlfli 


TiUl/i2 J Jl 
For |S<mi< ro, 0<Jna<|/3 
2 




^( 7 ->?.)+? (I- fx) 

/<a 


(64c) 

’ll) 

(84d) 




r ^u(f, ’?)— 


V;8*— m 2 * 

’T/11/X2 J Jr 


^(’l-’ll)+^(f-«l) 

ro* 


Ki^uVi) 


(64e) 


2(»ii-f:m2) 




TT/illdl 

For P<mi< 00 , /3<m2< “ 


^ (’i - ’ll) (?-?.) 

I Kiiu ’ll) 


l’o° 


(64f) 


C, 


2(mi+m2) 

vniidi 


If 


^(’i-’ii)+^(?-S0 
Ml 




_M2 
3 


’ll)«^flrfll 
(64^ 


The elementary-vortex distribution. — The reduction of 
equation (52), the formula expressing the potential due to a 
sheet of elementary horseshoe vortices, to the plane of the 
wing is of little interest since if Au is known »>a_o can be 


obtained by the simple relation <pt^o- 



Au dx. 


However, 


if the derivative with respect to 2 is determined on both 
sides of equation (52) and the limit is found as 2 approaches 
zero, the vertical induced velocity in the plane of the wing 
will be related to the vortex strengtli tliere. In a more 
physical sense, this will relate the slope of the lifting surface 
to the load distribution it supports. Such an equation is of 
basic importance. 

The mathematical expression to be studied is 


w=liin — CC^ Ap(^i,’ii) 

bz ivUiMi JJ 2 

T 


d^diji 


(65) 


The evaluation of w can be divided into two steps; First, 
the procedure necessary in order to carry the derivalive 
through the first integral and, second, the calculation of I 
where 


r-Htn + fA f— 

V-O 4tmiM2 j ^2 Jr rc ff 


d^drn 


(GO) 


Again the order of integration is important. As in the 
preceding discussion of the doublet sheet (equations (63) 
through (60)) and the source sheet (equations (61) through 
(64)), assume first that the area of integration is the one 
given in figure 2. Thus, mi is less than jS and theiji integra- 
tion is performed first. 

If e is introduced in order to evaluate the indeterminate 
form, equation (65) is expressible in the form 


«)=lini 

»-*0 

*-~»0 


(rni-l-mi^Fo b zy^ Ap 

47TMIM2 j-» ^'jxt * Te 2 


(67) 


where the teims Xi and X 2 are the roots of the qiradratic 
re*=0 and where the limiting process for e must be perfoi’ined 
before that for 2 . As before, the prime^e^ression denotes 
values for the particular case and eciuatiou 

(67) reduces to 


w=hm 

*-i0 

-.-.0 



[mi+m^Vu pa' £r/ 
4 ’TMiM3 



drii-\-I 


(08) 


By means of the mean value theorem, the first term on the 
right-hand side of the equation can be simplified. The pro- 
cedure involved in such an analysis was outlined in tin* 
derivation of equation (57). The process used here is iden- 
tical and equation (68) becomes 


T-’oVi8*-mi*Ap(^,7)) . J 

4 ^ 


( 69 ) 


The term /, defined by equation (66), can be expr«*sscd in 
a simple form by introducing the notation for the general- 
ized principal part (see equation (27)). This temi becomes 


,^ (mi + mj)Fo j:^ rp Ap(^i,i;d 


(70) 


If the integration of the above expressions had been taken 
in the opposite order or if the ra?ige of nii and f/ij had lieeu 
different, the residual term would change. This phenomena 
has been presented in connection with both doublet and 
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source distributions and is by now familiar. Finally, there- 
fore, expressions relating the slope of the vring surface to its 
loading can be written: 

For 0 <mi<^, 0<ma<j3 


To' 

4 2 

4iTllifl2 J 


V) 3 *— m2* Ap(i 77) 

t mi-f-miX 


4 2 

4 vjUl 7 t 2 J 

For 0 <mi<iS, 

CO 

1 C 

Vd*— nil* Ap(S, 17) 

miH-msJl 

Fo' 

4 2 

4 lT 71 l/l 2 J 

= 

j Jt 

j'oAp(^, 171) 
ro 2 

For [3 <mi< 00, 0 <m 2<|3 

w _ 

Vd*— Wi* Ap(f, 17) 

mid-miJl 


(71a) 


.7 jv *'oA2)(&,i}0 


Fib) 

^?l) 

(71c) 

(71d) 




Vq Ag(^i ^7i) 

ro 1 


(71e) 

(71f) 


J Jr fo 2 
For /3<wii< 03, iS<m 2 < 00 

Vq 4T/Ii/i2jrJ ro 2 

The doublet distribution. — ^The pertinent problem in this 
case is to finH the vertical induced velocity in the plane of 
the -wing as a function of the jump in potential across the 
plane. 

Consider equations (60) and tate the partial derivative of 
both sides with respect to 2 ; then find the limit of the resulting 
expression as 2 approaches zero. If equation (60a) is used, 
for example, there results for the first term 

y 5 1.,,. ,,1 , 1 /5A¥>\ 

5 J 2 A, (e., , j- j - 5 J +2 -gj 

which becomes 
—1 


(72) 


2 (mi-(-mj)Vi3-— m-i* 

^juaOS*— mi*) i7)+fii(mimi-[-d*) ^Ap(f,^)J 

The second term can be written 

i_o 5s 2vAtiAt2 X Jt re 

and this reduces to 


— jS*(mi+ma) C 

2viiiiii J 


r 

d^y^dm ^ 


A<g(Si,)?i) 

r 3 
^0 


2'n*Mi/i2 z-*o oz J Jr Tc 


Since, by equation (17), 

'X 2 


£ 


dtii 


[(X2-^0()?i-Xi)P' 


3 O 


the second term in expression (73) vanishes. Fin a lly, there- 
fore, the vertical induced velocity in the plane of the xring 
becomes: 

For 0 <mi<|S, 0 <m 2 <iS 

w=—— 1 — V 

2 (mi -bmi) V]S*— mi* L 5 f 

III (mimj-l- j3^ ^ Atf> (f, ij) J— 


Ap(fc, ^ 1 ) 




2tiii112 

^ , , - r Ms(mim34-i3*)^A¥j(f, 1 ;) + . 

)V/3— m 2 *L Of 


(74a) 


2(mi-hmj) 

;il(j3*— m 3 *) ^ A« 5 (f, n)J- 
J"or 0 < mi<^, ^ <m,i< « 


d^(mi+m3) 

2a-MiMs 




A¥>(fi, 17 0 


ro“ 


(74b) 


w== 




2(mi-i-77l2) 

;ti(mimj+i3*)^A¥s(f, 17 ) J— 


A«3(fe,7?l) 


2xj^iii2 


^TTHlHi J Jt 


ro“ 

Ag> (^, ^7i) 

ro* 


(74c) 

(74d) 


For ^<mi< “>,0<m2<[j3 
1 


■U) = 


2 (mi -[-m 2 ) 

;il(jS*— m2*)^A¥>(f, lj)J— 


fpiCminij+iS*) ^As>C|, ii)+ 


A?<&, 171) 


2vftlft2 


jg*(mi-[-m2) 


ro* 


2rilill2 

For ]3<mi< «»,)3 “ 




A?<&, 171) 


ro* 


d*(mi-f?n 2 ) 1 

Cj 

r A«<?i, 17 O 

2 vjLtlM 2 II 

g 

1 ro* 


(74e) 
(74f) 

d^dvi (74g) 
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SUMMARY OF RESULTS FOR r. C„ AND, u> IN THE PLANE OF THE WING; 
CARTESIAN-COORDINATE SYSTEM 

The special forms of equations (64), (71), and (74) when 
the J,i 7 axes become the Cartesian axes are given in the 
following sections. In these cases, ri~^y> wii— >0, and 

J7l2-> 00 . 


Potential in terms of vertical velocity, nonlifting case. — 


\ f dx, dy, (75) 


Pressure coefficient in terms of surface slope, nonlifting 
case. — 


O =— X ix v)——-^dxi-^dv-i a^i)^u(gi, l/i) 


(76a) 


Vertical velocity in terms of loading coefficient, lifting 
case. — 

v\ /3 Ap(g, y) 

Fo“ 4"1 ^ 


Ay(gi, y\) 


IX, X, <1 ■ 

{K-vO'V(l-i,)'-^>(v-vO’ 


X 


(y— yi)’V(a:— »i)*— )3%— 2/i)* 


(77a) 


"(77b) 


Vertical velocity in terms of surface potential, lifting 
case. — 


iQ /* 

Wu=—^ Au{x,y)—^;4-dxj-j-dyi 


<P»(giil/i) 

[(x — x,Y—P^(y—yi)^*'^ 


(78a) 



Vujx i, yi) 

[(x— r,)*— /3*(y— yOT^^ 


(78b) 


SUMMARY OF RESULTS FOR C,^ AND w IN THE PLANE OF THE WING; 

CHARACTERISTIC COORDINATE SYSTEM 

The special fonns of equations (64), (71), and (74) when 
the f, f! axes become the characteristic axes are given in 
the following sections. In these cases, f->r, 
and 

Potential in terms of vertical velocity, nonlifting case — 


<P 


1 f r Wu(ri,gi) 
vMoJtJ y(r— r,)(s— s,) 


dndsi 


(79) 


Pressure coefficient in terms of surface slope, nonlifting 
case, — 


2tP 




Vertical velocity in terms of loading coefficient, lifting 
case. — — 


<3 XX (r-ri)+(a-8i) Ap(ri,«i) 

Vo 4irXX [(«— Si) — (?•— ^’i)]’V('’— ?'i)(s — «i) 2 * * 

(81) 


Vertical velocity in terms of surface potential, lifting 
case. — 

■^^0 r r j j /'co\ 

4tXJ [(r-r,)(s-si)r ^ ^ 

PART II— THE DIRECT PROBLEM 


DISCUSSION 

The term “direct problem” sliaU be defined herein as a prob- 
lem requiring for its solution the evaluation of integrals with 
known integrands. Tlie three fundamental formulas pre- 
sented in part I, equations (64), (71), and (74), apply, respeiv 
tively, to source, vortex, and doublet distributions, and a 
consideration of them shows tliat there arc essentially only 
two different boundary-value problems of wing theory that 
lead to the direct classification. In the first of these problems 
(see equation (64)), the pressure coefficient is given by an 
integral involving the sliape of a wing having thickness, but 
no angle of attack, twist, or camber. The other direct prob- 
lem is represented by equations (71) and (74), where the 
angle of attack, twist, or camber of a wing having no thickness 
is given in terms of an integration involving, respectively', the 
wing loading or its streamwise integral, the discontinuity in 
-velocity potential. Tlie circumstances of tlie particular prob- 
lem will determine whicli of tlie tw'o alternative formula.*! is 
to be used. 

THE AERODYNAMIC PROBLEM 

The statement of the tw'o problems can be given from a 
physical vie\ypoint as follow’s. 

The thickness case. — The thickness of a wing that is sym- 
metrical above and below a horizontal plane is given and the 
pressure distribution over the wing is to be determined. Such 
problems are of special interest in the study of w’ings in a 
supersonic flow’ since tlieir evaluation is necessary for Uic 
calculation of the wave drag. 

The lifting case, — The load distribution on a lifting plane, 
a surface without tliickness, is given and tlie slope of the 
surface that w’Ul support such a loading is to be determined. 

THE MATHEMATICAL PROBLEM 

The mathematical statement of the two prolilems can bo 
made by referring to the equations exjiressing the three fun- 
damental formulas in the plane of the wing. For the thick- 
ness case, equations (41) or (64) apply, X„ is given over the 
wing plan form and is to be determined. For the lifting 
case, equations (43), (44), (71) or (74) apply, Ap/j or A^ is 
given over the wing plan form and Wu is to be deteimincd. 

The solution of problems in the direct classification depends 
only on the analyst’s facility in evaluating integrals. Although 
the integrations may be quite difficult to perform, this never- 
theless must be regarded as a question of teclmique and, in 
a mathematical sense, a direct problem is solved. 
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PAKT ni— THE INVERSE PROBLEM 
INTRODUCTION 

The term “inverse problem” shall be defined herein as a 
problem requiring for its solution the inversion of an integral 
equ^Ltion. In application to the study of problems in aero- 
dynamic wing theory, two different boundary-value prob- 
lems appear in the inverse classification. These are provided, 
as was discussed in the presentation of the direct problem, 
by' the two basic relationships that exist in the three funda- 
mental formulas. In those equations the two basic relation- 
ships are: First, the pressure is given in terms of an integra- 
tion involving the shape of a wing having thickness but with 
no angle of attack, twist, or camber; second, the angle of 
attack, twist, and camber of a wing having no thickness is 
given in terms of an integration involving either the wing 
loading or the discontinuity in the velocity potential. 

THE AEEODYNA2HIC PBOBLEM 

The physical interpretation of the two types of inverse 
boundary-value problems is made as follows. 

The thickness case. — The pressure distribution over a 
wing that is symmetrical above and below a horizontal plane 
is given and the shape of the wing is to be determined. To 
this bare statement of the problem, however, must be added 
certain auxiliary considerations. For example, it is physi- 
cally evident that solutions yielding wings with negative 
volumes must be excluded. Consideration must also be 
given to the question of wing closure. It is apparent that 
these tw'o conditions will serve to restrict the arbitrariness 
of the pressure distributions which can be prescribed. Fi- 
nally, the question of the uniqueness of the wing shape arises. 
For example, it is known that the thin-airf oil-theory solution 
in the two-dimensional case is imique, provided the prescribed 
pressure distribution is one leading to a real and closed wing 
section. In the supersonic three-dimensional case, however, 
these conditions are no longer sufficient to guarantee a unique 
shape from a given pressure distribution (although the re- 
verse is always true, L e., a given shape produces a imique 
pressure). This fact will be illustrated later (Part IV) in 
connection with quasi-conical flow problems. 

The lifting case. — ^The slope of, a lifting plate, a surface 
without thickness, is given and the resulting load distribution 
is to be determined. To insure imiquen'ess in problems of 
this type it is sometimes necessary to impose an additional 
condition. For example, it is necessary to assume that the 
Kutta condition applies to all trailing edges for which the 
normal component of the free-stream velocity is subsonic. 

THE IHATHEMATICU. PBOBLEM 

The mathematical statement of the two problems can be 
made at once. Thus, for the thickness case, equations (41) 
or (64) apply, Cj, is given over the area occupied by the wing 
plan form and Xa is to be determined. For the lifting case, 
equations (43), (44), (71), or (74) apply, ■«?„ is given over the 
wing plan form and Ap/? is to be determined. 

Of comrse, by definition, the solutions of both of these 
problems require the inversion of an integral equation. 
Further, these particular equations are known as singular 
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iut^al equations. Complete inversions to all the cases 
considered have not as yet been obtained. Some progress 
has been made, however, and the following section outlines 
one method by means of which certain singular integral 
equations can be inverted. 

ON THE INVERSION OF SINGULAR INTEGRAL EQUATIONS 

DEEnrmoNS 

An integral equation. — Consider the equation 

g(x) K(x, y)dx=w{y) (83) 

If v}{y) and iP(a;, y) are given functions and g{x) is unknown, 
equation (83) is known as integral equation, and more 
specifically as an integral equation of the first kind. The 
path of integration Li lies along the x axis (in this report only 
real variables are considered although the methods and 
results can be generalized to include complex variables) and, 
in general, can depend on y. The term K(^, y) is known as 
the kernel of the integral equation. 

A singular integral equation. — ^An integral equation is 
referred to as singular either when the path of integration, Li, 
has infinite extent, or when the kernel, K(x, y), is infinite at 
points of the interval ii. In other words, equation (83) is a 
singular integral equation if KQe, y) is unbounded somewhere 
onii. 

An integral transform. — .Again consider equation (83). If 
both sides of this equation are multiplied by the function 
jff(X, y) and integrated with respect to y along the interval 
Li (which is, of course, independent of y but can be a function 
of X), the equation is said to have been transformed and 
the operator 

y)dy (84) 

is referred to as an integral transform. The resulting 
egression 

dxgix) H(\, y) K(x, y)= j^wiy)Hi\ y)dy (85) 

is obviously a function only of X, both x and y being dummy 
variables of integration. 

mHEBENT SINGCLAHmES 

An inherent singularity can be defined first in terms of a 
function of two variables. Consider the fimction/(x, y) and 
let the point a, h lie somewhere in the x, y plane. Then an 
inherent singularity will be said to exist at the point a, b if 

lim e*/(e|-[-ffi, eij-j-J) 

e — *0 

where 

x—a=e^, y—b=et] 

In other words, a square of width 2e is first placed on the x, y 
plane with a, b at its center, the function/(x, y) is then eval- 
uated at any point on the boundary of the square and multi- 
plied by one-fourth the area of the square. Finally, in the 
limit as the width of the square vanishes, if this product is not 
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zero, the function /(at, y) contains an inherent singularity at 
the point a, 6. ... 

Such a concept can obviously be_ generalized to include 
functions of three and more variables. For example, the 
function /(?!, Xi, x^) contains an inherent singularity at the 
point di, Oj, a„if 

lim + , e?a+cf«)?^0 

e-O 

where 

x,—at=6^t 

RESIDUALS 


Consider a double integration with respect to x and y of the 
function /(a:,y) over the area S in the x, y plane. Perform the 
integration of the same function over the same area but with 
the order of integrations reversed. The difference between 
the results of these two operations will be defined as the 
residual. Thus 

^dy^dx fix, y)— ^dx^dyj{x,y) (86) 

Ordinaiily the residual i? is zero, since the order of inte- 
gi-ation for a double integral is usually immaterial. In the 
manipulation of singular integi-als and singular integral 
transforms, however, a nonvanishing residual often exists. 
The evaluation of the residual can be accomplished in the 
following manner. Let the point di, hi be an inherent singu- 
larity in the area S. Then the residual from such a point is 


Setting ef-f-di=a: and one can wite 



df lim eri + hi) 


dft ]imt^f(€^+ai,eri + hi) 


(87) 


Hence the necessary condition for the e.xistence of a residual 
is the occurrence of an inherent singularity in the area of 
integi’ation over which the double integration is performed. 
The total residual is the sum of the residuals from each inlier- 
ent singularity in tlie area involved. 

THE NULL TRANSFORM 

Definition. — The integi-al operator H{\, y) dy is said 

to be a null transform of order n to the function K{x, y) in 
the interval Z* if 


where x, X and, of course, y are on Zj. 
Examples. — The operator 




dy 

i\-yr^ ■ 


( 88 ) 


is a null transform of order zero to the function 


K{x,y)=lK^z 


in the interval r<y<X. Thus, (sec equation (18)), 

-p-j=£M==0 ' (89) 

\(y-x) ix-yr 


The operator 


[■ m\,y)dy= 


X) 


is a null transform of order zero to the function K(x, y) = 
ll(y—x) in the interval X<y<l. Thus 

(y-3-)V(i-y)(x-y) 

Finally, it can be shown from equation (90) that 

where a and b are constants, is a null transform of order one 
to the function K(x, y) = l/(y—x) in tlu'. interval a<Cy<ih. 
Thus 


TUB INVERSION OF SINGULAR INTEGRAL EQUATIONS BY 
MEANS or NULL TRANSPOR.MS 

Consider an integral equation of the fimt kind 

‘^iy)=Jj^gix)Kix,y)dx (92) 

such that the kernel K{x, y) tends to infinity’ as x n]>proachcs 
y and let the point x—y lie in L\. Equation (92) is, by 
definition, a singular integral equation of the first kind. 
Apply to both sides of this equation the integral Iransforni 


(y-x){\—y) " (d 


d<ar<6l 

<\<b) 


(91) 


J'^H(X, 2 /)dy so that 

J^w(y)fZ(X, y)dy— dx^(3r)//(X, y)K(x, y) (93) 


Suppose that the area ZiH-Z* of the double hilegral is 
bounded by a simple closed curve having the property that 
any line parallel to the x or y axis crosses its boundary at 
most twice. For such an area, it is alwaj's possible to write 
the revemed form of the double integral in equation (93) as 

^^f^^j^dygix)Ui\y)K{x,y) (94) 

w’here Lt can be a function of x and X, and Li can be a function 
only of X. 

Subtracting expression (94) from the double integral in 
equation (93), one finds 
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y)K{x, y)= 

\^dx^^iy g{x)H(X, y)K{x., y) +i? (\) (95) 

where R{\) is, by definitioa, the residual. Hence, equation 
(93) can be rewritten in the form 

J^!£?(jf)fl'(X, y)dy= 

f dx r dyg{x)H{\. ymx, y)+B{\) (96) 

J -^4 J ^ 

One can now show that if y)K{x, y) contains an 
inherent siugularitj' at the point x=y=\, equation (96) is 
the inversion to equation (92) when y) is a null trar^- 
foim of order one or zero to the kernel K (a:, y) in the interval 

First, it is necessary to relate 1?(X) to ^(X). By the denm- 
tion given as equation (87), .ff(X) can be written 

i?fX)=ffo(XH“ [ dj] f dfZ‘im^e$-{-X)iiL(e5+X,ei;+X)fl(X,«5;+X) — 
J — 1 J -1 «-*o 

r cff r di7fim.^(e^+X)liL(«f+X, ei;+X)H(X, ei;-{-X) (97) 
J-l J-l e-K) 

where i?o(X) is the sum of the remaining residuals (if there 
are any) from the other inherent singularities that might 
exist in the area Li-\-L 2 . This reduces to 

R(\)=E,i\)+g(\)R*i\) (98) 

where 

i?*(X)= { dti{ df /tmif(e^+X, e7j-t-X)H(X, ei7+X) — 

J-l J-l ^ 

(* df r dijZtm i!L(e^-i-X, ei7-f-X)H(X, 6 t;-|-X) (99) 

J-l J-l €-» 

If inherent singularities other than the one at x—y—\ exist 
in ii+ij, they must be at a point on the line x=y. The 
residual at such a point say x=y=a, would be the product 
of g (a) and J?o*(tt)> the differepce between the two appro- 
priate double integrals. Hence, i?o(X) cannot contain g(k). 

Now, if y) dy is a null transform of zero order, 

equation (96) becomes 

f w{y)m\ y)dy=Roi\)+g(\)R*W 

.1 ^2 

or 

5(x)=-gJ^ 

which is an inversion of equation (92). Further, if 

^Jd(^,y)dy 

is a null ti’ansform of the first order, so that 

[* H(X, y) K(x, y)dy=C(\y 
J -^3 

equation (96) becomes 
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J^iij(i/)H(X; y)dy=C(\)j^gix)dx+Ro(\)+g(\) R*{\) (101) 

If Li is independent of X this already is an inversion of equa- 
tion (92). However, if Li contains X then, after rlividing 
through by C(\) and . taking the derivative with respect to 
X, a first-order differential equation in g(\) results. This is 
considered to be an inversion. 


THE INVERSION OF SOME PARTICULAR SINGULAR 
INTEGRAL EQUATIONS 

ABEL’S INTEGRAL EQUATION 

Consider the special form of Abel’s integral equation 


rim. «<» 

Ja \y — X 

It has been shown (see equation (89)) that 


( 102 ) 


f*’ dy 

X (X-2/)*'* 


is a null transform of oiffer zero to the kernal l/vy — ^ Ih 
the interval x<y<\. Hence, applying the transfoi-m 


r*- dy 

to both sides of equation (102), revei-sing the order of inte- 
gration, and noting that 


-f dy f 

Jtt. Ja 


dx 


iPd) 


-^y—x (X— y)»/- • 

■y g(x)dx-^ dy , ^ — -t-l?(X)=i?(X) 

J. Jx yV—x (X— y)®/* 

(103) 

leads one to the result 


The only inherent singularity that appears in the area of 
integration occurs at the point x=y=\. The residual is 
obtained by integrating over the shadetl area in figure 5 and 
finding the Hmit as e goes to zero. 



Figcse 6.— Eegtoa oflntceratlon Cor equation ClOJa). 
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Thus 


l?(X)=limr-^ dyf dx 

yy—^O^—yr 

r g{x)dx-f dy ■ ■ J (lOSa) 

Jx-« ,h y'y—x{\—yy‘*j 

and, since the second double integral is zero, the transfor- 
mations ef=X— a: and eri=\—y reduce this to 

B (X)=Iim -f dn rVf 

t-io ji ji Y ^ — n 17*'“ 

Finally, in the limit as e goes to zero, 

2Ti(X) (105b) 

Substitute equation (106b) into equation (104) and the 
inyersion of the integral equation (102) can be written 




(106) 


By applying the definition of the finite part, one can rewrite 
equation (106) in the alternative form 

which is the form of the inversion usually presented. . 

THE AIRFOIL EQUATION 

The study of the singular integral equation known as the 
airfoil equation is closely associated with the study of bound- 
ary-value problems related to Laplace's equation in two 
dimensions. These boundary conditions are sometimes 
given along a straight line as is the case, for example, in the 
linearized study of two-dimensional subsonic wings. If the 
boundary conditions are given along a suitably prescribed 
curve, the curve can, by application of complex variable 
concepts, be mapped onto a straight line. For example, the 
JoukowsM transformation maps a circle in one plane onto 
a straight hne in another and in both planes the governing 
formula is the two-dimensional form of Laplace’s equation. 
The solution to such boundary-value problems can be reduced 
to the inversion of the following singular integral equation: 


w(y)-. 


p g(x)dx 
Ja y — 


a<y<b 


(108) 


where a and b are constants. 

It has been shown (see equation (91)) that 


X 


’’ yKb-yXy-a) , 
\-y ^ 


is a null transform of order one to the kernel Ijiy—x) in the 
interval a<y<b. Applying this transform to equation 
(108), and using the definition of the residual, one obtains 


p «?(y)V(&-i/)(j/-g) 




which reduces to the form 


r‘ wiiiU(b-y)(v-c ) j r* 1,09) 

Ja A— y Ja 

Again the only inherent singularity in the area of in- 
tegration is at the point x=y=\. Evaluating the residual 
according to equation (87), one finds 

s(x)=«(x)v(j-x)(x-a)[/;</,j;;^- 

C - 
J-t 


and since 


dv _-l 1-e 

■'”T+£ 


this becomes 


J?(X)=2^(X) V(6-X) 




( 110 ) 


By the combination of equation (110) with (109), the inver- 
sion to the airfoil integral equation (108) thus becomes 

^(^= : — 


(Ill) 

It is apparent that the inversion to equation (108) provided 
by equation (111) is not unique (because of the existence of 

the term- J* g{x)dx which can be thought of as an arbitrary 

constant). Hence, in the application of equation (108) to 
physical boundary-value problems it is not sufficient to 
specify the value of V) along the y axis; some additional con- 
dition must also be supplied. E.xamples of such additional 
conditions in the study of aerodynamic problems are the 
specification of closure in the study of two-dimensional sec- 
tions and the assumption of the Kutta condition along the 
trailing edge of two-dimensional lifting surfaces. 

THE SUPERSONIC DOUBLET EQUATION 

The general concepts of the method just applied to the 
solution of single-integral equations with a singular kernel 
can also be used to invert double-integral equations with 
singular kernels. Success in solving these more complicated 
forms depends again on the discovery of an appropriate 
integral transform — in this case a double integral transform— 
and the usefulness of these operators depends, in turn, on 
both the structure of their integrand and, what is just ar 
important, the space (now four-dimensional) of integration. 
As it turns out, however, inversions can be obtained in many 
cases that are of importance in the study of supersonic 
aerodynamic problems. 

Case 1 — Supersonic leading edge. — Consider the equation 
that gives the vertical induced velocity in the plane of the 
wing in terms of the perturbation velocity potential on the 
upper surface of a lifting surface in a supersonic free stream, 
equation (82), 
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W„(f2,S2) 


-Mo 1 

pi 

f «5«(ri,si)cfridsi 

4a- J 

n 

r [(ra— ri)(s2— Si)]®/* 


Tin, «j) 


( 112 ) 


JJ 

rocr, 0 


drjdSj _ 2j9 
iJir—r^(s—s^ ^0 


// 

raff, a) 


dvj dsj 
ro(r,s]n,sii 


Cl 13) 


The area r(r 2 ,fi 2 ) represents the area ■within the forecone 
from the point Ti,S 2 at which the vertical induced velocity 
is being measured. The details of the solution to equation 
(112) will now be presented for two different types of bound- 
aries to the area r(?* 2 ,S 2 )) is, from an aerodynamic 
standpoint, for two different types of ■wing plan forma. 

First consider the case when r=ro is an area such as the 
one shown in figure 6. The two lines ri=r 2 and Si =« 2 , which 
represent the traces of the hlach forecones from the point 
Tt, Si, form two bounds of the area while the third is given 
by the ■wing leading edge, the equation for which may be 
written as either Si=/oC?*i) or J*i=/o*(»i)- 


Vo 




(a) n. III plane. 

0>) n, «i plane. 

FinuBE 6.— Segfon of btegiatfon fpr wbg wlUi snposonic leading edge. 


In the present case, it is assumed, furthermore, that «i=/o(ri) 
is a continuous monotonic function ■with a negative slope, or, 
physically, that the ■wing has a supersonic type leading 
edge.® 

The kernel of equation (112) is formed by the product of 
two functions, one independent of ri and ?*2 and the other 
independent of Si and S 2 - It follows, therefore, that if an 
integral operator is used that contains the product of two 
linear functions of «2 and r 2 to the —1/2 power, there exists 
the possibility of obtaining a null transform of order zero 
according to the equality given as equation (18). One can 
.show, in fact, that the operator 

‘ EapeisonCo and sabsonla edges have the propert:r that the normal component of the tcee- 
stieam velocltr S Enpets(mla and snbscHiIc, cespe<^Tel7. 


is a null transform of order zero to the right-hand side of 
equation (112) where Ta(r, s) is the same area that appears 
in the integral equation. 

If the transform given by equation (113) is applied to 
equation (112) the resrdting expression is 





WudrjdSi 

I'd 


irMo* J J J J ro(r,s,-r2,s*)ro®(r2,fi2;ri,Si) 

ri(r,a) rt(ri »i) 


(114) 


Consider nest the ri^t-hand side of this equation but 
with the r 2 , Sg integrals taken first. Then by d efini tion of 
the residual 


nr Mo* J J J J ro(r, s ; r^, Sg)ro*(r2, ri, sd 

Ti(r, ») T,(n, J,) 


— 2;8® 
vHfo*. 

R(r,s) 


fj 


VuiruS^dridsi 
[n(r, »)+ro(n. »j)T 


ff 


drjdSi - ■ 

ro(r,s; Ti, aa)'*o*('’i, SzJ ^i, 

(116) 


When the integration is made first with respect to and 
« 2 , the area of integration for these two variables can be 
■visualized with the aid of figure 7. In such a case the 
points ri, Si and r, s are fixed and the area is simply the 
one which lies in the forecone from the point r, s and in 
the aftercone from the point ri, Si. This is represented 
by the shaded area in the figure. It is apparent from a 
study of figure 7 that, when the edge S 2 =/o(r 2 ) is a mono- 
tonic function as pre^viously defined, the and s% integrals 
are always taken between the limits ri and r and Si and s, 
respectively. Hence, according to equation (18) the inte- 



Fiquee T.— Area of integration for fired n, at, «; equation (116). 



1286 


REPORT 1054 — ^NATIONAL ADVISORY COMMITTEE FOR AERDNAXJTICS 


*p'al term on the right-hand side of equation (115) vanishes 
and equation (114) becomes 

J |’^^^|^=7?(r,s) (116) 

To(r, ») 

The complete inversion of the integral equation can now 
be realized if an expression for R(r,s) can be obtained in 
terms of ^«(r, s). The evaluation of this residual term 
follows along lines quite similar to those used in calculat- 
ing the residual for the special form of Abel’s integral equa- 
tion given previously. Now, however, there is no longer 
a single inherent singularity; rather, the lines ri=rt=r and 
are densely covered with them. First consider 
an integration made over the region close to the line «i= 
8i—fi (i. e., the sum of the areas a and b shown in fig. 8). 



Figure 8.— Region of integration for evaluation of residual, equation (116). 


Make the substitutions Si=s— e<ri and sj=s— etra and take 
the limit as e goes to zero. There results for the portion 
of the residual due to the inherent singularities Ra+tir, d) 
along the line .9i=S2=s the equation * 
f^a^6(^,s) = 


8irfi 





dffi 


s) 


(117) 


'I’his becomes, after integrating with respect to vi and o-j 
(which can bo done immediately since the limits of the a 
integral do not now contain o-j and, further, since the r,, 
plane contains no inlierent singularity) 

(The drat term on the right-hand aide of equation (US) oontrlbutca nothing to the residual 
since it always contains an Integral equivalent to that in equation (18) and, therefore. Is 
Identically tero. This same phenomenon also appeared In the study of the residual appear- 
ing In the Inversion of Abel’s Integral equation. (See the development of equation (105).) 


R.Ur,s)=^ f dr, 

4)3 J/o‘(» *J/o*(.) (rj- 


Vu(.ri, >i) 




( 118 ; 


It is now proposed to reverse the order of integration of 
the double integral in equation (118). But the r,, plane 
contahis an inherent singularity at the point rj=ri=r. By 
definition of the residual fram this singularity as R'(,r,8), it 
follows that 




R'[r,s)=R'{r,s) 


since the integral term is zero by the equality that has been 
used repeatedly. The evaluation of R'{r,s) follows the iden- 
tical line of argument used in obtaining equation (105b) from 
(105a). Hence, setting ra=r— and ri=r— ep,, and letting 
ego to zero, one finds 


Ra+i,(r,s)=R'ir,s)= 



<Pu{r, s) 



dpi 

{Pi— 


2)3 


vuir, 8) 
(119) 


The part of the residual in equation (IIG) contriiiuted by 
the inlierent singularities in the areas a and h in figure 8 is, 

therefore, — «o«(r , »). A similar calculation shows that the 


singularities in the areas c-+ 6 give the same result; and, finally, 
a calculation for the ar(>a 6 itself also yields — p«(C) f). 

Tlie value of 7? inequation (1 10) isobtaiiiid by comliining 
the results for the various arc'as. Hence, 


R=-Ra+^-\-Rl>+c — Ri 


-vAV 

2^ 


<Pu{r, s) 


( 120 ) 


Combining the result e.xprcssed by the last equation witli 
equation (116), one can finally write for (he inversion of the 
integral equation (112), when t=Th the e.xpressiou 


<Pu{r, s) 


— 1 ('I' Wudrjdst 


( 121 ) 


Case 2 — Combined subsonic and supersonic leading 
edges. — Consider the equation 


— A J o r C v>v(ti-, Si)Jri</S i 

4ir 77/ [(^2— fi)(Sj— Si)]‘^* 

Ti(rj,«) 


( 122 ) 


where ri(r 2 , Sa) is the more complicated area shown in figure 9. 
Again the lines ri=r 2 and 5‘i=Sa, which are the traces of the 
forecone from the point Si, form tAvo bounds of the area. 
The remaining two boundaries arc formed by the curves 
8i==/„(/’i) (or in the inverse sense ri=Jo*{si)) and Si=/i(ri) 
(or ri=/i*(si)) where /» has the same definition it had in llie 
study of Case 1, that is, a monotonic curve with a negative 
slope. The curve Si=/i(ri) is also a monotonic function, but 
with a positive slope. For convenience the origin is placed 
at the point of intersection of the/, and /, curves. 
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(a) A, III plane. 

(b) fri, Ji plane. 

Fiocbe 9. — Beglon d IntegncUon tor Ting Tritk combined subsonic and supersonic edges. 

A physical interpretatioDi of the area t, with regard to prob- 
lems in supersonic wing theory is simple. The lines ri=ri 
and Si=S 2 i as has been mentioned, are the traces of the Mach 
forecone from the point rj, Sj. The/o and/i curves are the 
edges of the wing plan form, the line Si=/o’(ri) representing a 
supersonic leading edge and the line ■Si=/i(vi) representing a 
subsonic leading edge. 

Proceeding exactly as in Case 1, consider the operator 
f f drids2 

and find whether or not it is a null transform of zero order to 
the kernel of equation (122), that is, whether or not the 
equality 


0' 4»«(''i, s^dridsij^ 


rtCt,rl+Ti(n.n)1 



(123) 


is satisfied. A studj' of figure (10) shows that it is not, since 
for a point ri, Si located above the line ri=/i* (s), a portion of 
the Si integral becomes 


Jn (Si—Siy'\s—Si 



I 


J’s 

noxTBE 10.— Area ot Integcatioa for fixed ri, < 1 , r, i; eauatfoa (123) . 

which is not zoto unless /i(ri) is identically equal to s. Por 
the same reason, portions of the rg integral ^viIL not vanish for 
the point rj, Si so located. 

However, the construction of a null transform of zero order 
can be accomplished by studying the above failure. Thus, 
if the 8s integration were carried between the limits «i and s, 
for every location of the points ri, Si and r, s, the operator 
J'h(r, 8, rs)drgJ'is—Si) would be a null transform of zero 

order regardless of the choice of A. The area of integration 
producing such a transform is simplj' the one shown in figure 
11 (a), an area bounded by the lines S 2 =s, Ti=r, 8i=foirg) and 
rg=fi*(8) where ri=Js*(8) is any line such that/ 2 *(s) >/i*(a). 
This area shall be designated as Tg(r, s). When such a trans- 
form is applied and the S 2 integration performed first, it is 
apparent from figure 11(b) that the limits of the St integral 



(a) Area of tutegratfon for transform. 

Cb) Area of Integratiou for fixed ri, ii r, t. 
Fir.uBE IL— Integratiou regions for equations (120. 
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are always and s, even, when the point rj, Si is located above 
theIuierj=/:*C«)- Hence, if the residual i?(r, a) is defined by 
the relation 

Ma fw r dSi C C <Pu(TuS^dndsi , 

•' ’ n(n,«) 


’’Kr,*) 


[’lC'’.»)-H-iCr 2 , * 2)1 (124a) 

the second term on the left-hand side is zero, and equation 
(122) transforms into 

r h{r,s,TiidrS ^’‘^^ ■da.=R(r,a) (124b) 

J f2*W J/ofr2) Vs— Si 

The evaluation of the residual in equation (124b) follows 
the same pattern used in Case 1. Thus, after isolating the 
line of inherent singularities, R(r, s) is given by the expression 

5(r,s)= 

4:5tJ/2 


/2*W 


hdr. 


r* daj 
J »-« -yja — S 2 J /i’ 


da dr 

SiJ ^ ‘ jT-e ’ [(^2— r0(a2— »i)]®'* 


which becomes (after setting Si=s — « t 2 and Si=s — evi and 
letting e go to zero) 

R=^r hdr^-p^ (125) 

2 J/ 2 *W J/i*« (ri— ri)®" 

It is now a pparent that after setting / 2 *(s)=fi*(a) and 
h—lj^-r—ri, equation (125) can be evaluated exactly as 
equation (118) was evaluated, so that for these values of 
■fi* and h, R (r, s) = —Mq-k<pu (r, s) . Substituting into equation 
(124b), one finds for the inversion of equation (122) the result 






^ 2 ) 

rj) (s— Si) 


(126) 


where 7-3 is the shaded area in figure 12. 

Case S — Mixed boundary conditions. — ^Another very im- 
portant Jdnd of integral equation which can be solved directly 
by the proper choice of A in equation (125) is the “mixed” 
type problem, the boundary values of which are illustrated 
in figure 13 . In this particular problem is known over the 

portion- of the r, a plane bounded by the curves a=fa(r) and 
s=fi(r) (the curve a=f 2 (r) is a monotonic function with a 
positive slope just like »=/i(r)), while over another portion, 
bounded by the curves a=fi(r) and s=/i(r) the quantity 

-Uu(r,a)=^ ^ (This corresponds in 

aerodynamic applications to the specification of the vertical 
induced velocity over the former region and the loading over 
the latter.) It is also assumed that a) is continuous 
across the line a=ft(r) and that «j(r, a) vanishes along the 
line s=/i(r), that is, <p(Ji*ia), a)=0. 

Taking for the value of A the expression 


-vr 


— ra 


Va 



Vb 



Ca) ij,ifiplane. 

(b) fi,si plane. 

Fiocse ' 12.— Integration area r, (eauatbn (iseu. 






(aj z. If plane. 

G>) r, 8 plane. 

Fiocse 13.— Mixed boocdary conditions. 
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and for the area of integration the r^on n(r, s) shown in 
figure 11(a), one finds from equation (125) that the residual 
can be written 


2- J/a*C*) 


(ri, s) drt 
fi*W (r*— 


The partial derivative can be taken through the second 
integral sign (first integrating by parts in the case of 5 /Sa) 
since a)=0, so that 


B 


=Mo r 

2 J/a 


r* 


dri 




dr * 

/2*W ^J/r*C«) ^ fjr — — ri)®^* 


Part of the boundary condition is that Uu(r, s) is given 
over the region indicated in figure 13. Hence, the residual 
can be written 








Vr— ra(r2— ri)*/* 
drg -f * dri . .. f) 


i- 


The latter of these two terms is again just like the one given 
in equation (118) so its evaluation is immediate. The first 
term may be simplified by reversing the order of integration. 
(There is no residual since the point r=ri=r 2 is not included 
in the area of integration.) Finally, E becomes 


i2=2/3Vr-/2"(a)J^ 


•^**f*) ttg(ri, s)dri 


-2ffjS'Uu(r,s) 


(r— ri)V/a*(s)— ri 
and the inversion of this mixed type problem can be written 

v-JXu s)dr^ 


’ll Jr, s)-- 


TT 


X 




_i 
2)3 V 


(r— ri)V/2*(s)— J": 

r . drz / d a \ r* ' wjrz, Si)dSj 

J /2*(J) -v r— r,\3r2 ‘"da/ J /„( 


-y's — 5a 


(127) 


An equation which does not impose the condition that 
'Ua(r, s) be continuous across the line s=fjr) can readily be 
developed by using the operator h=\l-\jT — rj. This leads 
directly to the result 


f' 

vAfo Jf,* 


/ 2 *W (p„ (ri, s) dri 


A*(.) (r— ri) Jf*{s)—r^ 


0 f* is. 

J /ocn*} 


viu (?~2. sj 
’V(^— Sj) (?■— rj) 


(128) 


THE SCPEHSONIC SOUHCE EQUATION 

The null transforms to the supersonic doublet equation 
were constructed by appljrng an integral operator having 
an integrand which, when combined with the kernel, would 
produce the integral 

dx 




(r— 6) 

and having an area which, when combined with that of the 
integral equation and traversed in reverse order, would pro- 


duce the limits a and h on the integral. If one now considers 
the supersonic source equation (see equation (79)) 

.0-^' ff (129) 

- yA,V('*2— J*!) (si— «i) 


■rCrj,»d 


it is apparent that the null transform will differ from the 
transform used for the supersonic doublet equation only in 
the exponent of the intregrand. Thus, for the area ro(r 2 , Sj) 
(see fig. 6) the operator 


X 


I. 


dr.ds. 


_ J [(r— rs) (s— aj)p/* 

U{r,t) 


is a null transform of order zero to equation (129), and its 
application yields the solution 


wjr,s)-- 


Mq r r ¥>u(?*2, gg) dr^ Sj 
JJ [(?•— r^(s— s 


4ir 


[(r-r*)(s-s*)]='/* 

For the area rdr^. Sj) (see fig. 9), the operator 

dr ^ds 2 
J [(r— T 2 )(s— 4 


(130) 


X 

J 

’■iCr.*) 




is a null transform of order zero. Its application, under 
conditions like those specified in the development of equation 
(126), yields the inversion 




s^dvidst 


nix,*) 


/•=)(s-S2)P» 


(131) 


Equation (130) gives vertical velocity in terms of a prescribed 
surface potential for a wing with a supersonic leading edge, 
while equation (131) does tiie same for a wing with a leading 
edge which is partly subsonic and partly supersonic. The 
transform with the area vg (see fig. 12) can also be used to 
obtain inversions to equation (129) under conditions such as 
those imposed in the development of equations (127) 
and (128). 

DISCUSSION OF INVEBSION OF SUFEBSONIC SOURCE AND DOUBLET 
EQUATIONS FBOAI A PHYSICAL BASIS 

Each of the above examples used to filustrate the applica- 
tion of the null transform method to the inversion of double 
integral equations represents the solution of a class of super- 
sonic problems. However, most of these solutions are well 
known and were originallj^ obtained by reasoning that was 
suggested by knowledge of the physical structure of the 
problem. For example, the reciprocal relation between the 
source and doublet int^ral equations 

X Mo J r„ J y(j— rj) (s —S i) , 


,,, (r :zM®X r vjri,s;)dndsi 

4v Xo J [(r-ri)(s-sdr 




(132) 


when considered with respect to the area ro, has a simple 
physical meaning. By defibaition tq is an area bounded by 
the Mach forecone and a supersonic leading edge. From a 
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physical standpoint it is clear that the flow fleld at a point, 
affected only by a supersonic edge, on the upper surface of a 
wing cannot be influenced by the sliape of any pai-t of the 
lower surface of the wing. In (ffher words, the upper and 
lower surfaces are noninteracting. Hence, the upper surface 
of a wing does not “know” whether it is the upper surface of 
a lifting plate that is supporting loading and has no thickness, 
or the upper surface of a wing section that is symmetrical 
above and below, the 2=0 plane. Thus the source equation 
must be the inversion to the doublet equation and vice versa, 
and from a physical point of view the reciprocal relation 
given by equation (132) is obvious. 

The solution to wing problems involving one supersonic 
and one subsonic edge, giving an area of integration for the 
integral equations corresponding to the area u in the pre- 
vious discussion, was originally obtained by Eyvard (see 
reference 7). It should be noted here that the inversions to 
the source and doublet equations (equations (126) and (131)) 
considered with respect to the area rj, no longer foiin recip- 
rocal relations. 

Finally, the examples presented. herein.with regard to the 
mixed type of problem have also been deiived (see references 
8 and 9) using more or less physical arguments. The solu- 
tions to these mixed type problems form the basis of a lift 
cancellation technique that provides a very useful extension 
of Evvard’s original discovery. 

ITERATIVE METHODS OF SOLUTION 

Other types of plan forms. — The question that naturally 
arises from a practical viewpoint is how the source and 
doublet equations can be inverted when the area n is not of 
the two special kinds discussed, or, in otlier words, when the 
wing plan forms are complicated bj' having more than one 
monotonic (in the r, s plane) subsonic edge. The answer 
must be that, unless null transforms with respect to these 
new integration areas can be discovered, the methods dis- 
cussed here will not give the direct inversion to the problem. 
Several possibilities remain, however, so that even when the 
null transform cannot be found the concepts of the residual, 
inherent singularity, etc., can be used to simplify, if not 
solve, the supersonic source and doublet integral equations. 

With respect to the lift cancellation techniques already 
mentioned, references 8 and 9 outline these methods in 
considerable detail and show how they can be applied to find 
the loading on wings in regions affected by two or more 
subsonic edges. 

Regions influenced by multiple reflections of Mach 
waves. — A more direct example of how some of the concepts 
presented heretofore can be applied, even when the null 
transform is not available, is given by considering the fol- 
lowing problem: Find the loading at the point (r, y) on the 
flat wing tip shown in figure 14. If the r, s coordinate sys- 
tem is used and the operator 

r dra / d i ^ N p ds; _ 

J /a* (») J /a(r) 

(where the area of integration is region 4 in fig. 14 (b)) is 
applied to the doublet equation (equation (82)), two results 


Vo 




(b) n, h plane. 

FifiVBE 14.— Tip of swept forword wing and region* of Integritlon, 

can be anticipated: One, since the loading off the wing is 
zero, the residual will be— 2TjSWa(r,a) ; and the other, the trans- 
formed integral with the (in the notation of equation 
(124a)) integration performed first, is zero for all points ri,Ui 
lying in regions 2 and 4 in figure 14. These results follow 
directly from the discussion presented above in case 3 of 
the similar integral transform applied to the doublet equation 
(equation (122)) with the n area. Without proceeding 
fm-ther, therefore, it is apparent that the original integral 
equation has been reduced to: (1), an integral of the known 
functiou.ii;tt over region 4 in figure 14; and (2), an integral of 
the unknown function ipu over the regions 1 and 3 in figure 14. 
(It can be shown that the integration over region 3 will also 
vanish.) But regions 1 and 3 are aheail of the Mach forccone 
from (r,«). Hence, by repeating the above process for r^ions 
farther and farther up (toward the origin) the wing, the 
problem must eventually be reduced to one of finding the 
solution for an area such as ti; in other words, to a problem 
involving only two edges, one subsonic and one supersonic. 
Since the latter problem is solved, the one considered in this 
section i^Iso (theoretically at least) solved. 

Triangular plan form with subsonic leading edges. -As a 
flnal e.xemplification of the preceding concepts, consider the 
problem of finding the loading on a flat triangular wing tly- 
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ing at a supersonic speed but with both leading edges sub- 
sonic (see fig. 15). To the supersonic doublet equation 
(equation (82)) apply the integral transform 



where the area of integration is region 4 in figure 15(b). 
There results 

^ dr^ {' d8; _ 

J./m Jr/j. \(r—r^ (a— sd 

97. (ti, gi) dr^^dsi 

V(J"— ra) (a— So) [(r,— rO (a^— Si)l®^® 

(133) 






where Ti{r^, Sj) is the original area of integration as shown in 
figure 15(a). Since uv is a constant equal to —Vooc, the left- 
hand side of equation (133) becomes 




(b) 

(a) xt, t/i plane. 

(b) fj, »t plane. 

Eiocbe is.— I ntegration areas for trlinanlar \rfng. 


where ju=(l-[-mj8)/(l— mj8), m being the slope of the leading 
edge in the ar,y plane. Inverting the order of integration 
and evaluating the residual, one finds that equation (133) 
can be written 


9 uir, s)= 


4TV^(r-0(»-r) 




9u (Ti, 8 i ) dridsi 


(r-rO(a-ad^g-a,)(j-rx) 


where is the area shown as region 1 in figure 15(b). 

Equation (135) is not, of course, the inversion of the super- 
sonic doublet equation for the triangular flat plate. In fact, 
it simply represents the transformation of the doublet equa- 
tion, which is a singular integral equation of the first kind, 
to a singular integral equation of the second kind. How- 
ever, this latter form has the advantage that it is readilj' 
susceptible to the process of iteration. Thus, in the par- 
ticular case of equation (135), it is possible to take as a first 
approximation to the value 

and, as the second, the value 

[, 1 r r 1 


-MJ 


VF'O (»“'■') 


and so on. By means of the substitutions 

a-hfAtpi-l-ag-L *) 

^ ^ 1-f-pi-j-o-i I 

, r r (138) 

. rd-rpt-hSjXQ-i, ( 

^ ‘ l + Pi+vi J 

the double integral in the equation for 92 can be evaluated. 
Thus, after some manipulation, one finds ^ 


JJ 




=2x[2Ei— (1— IT® (139) 


where Ki and Ei are complete elliptic integrals of the first 
and second kinds, respectively, with moduli l*i=l/p. Using 
the identity 

r Notica that the area of Integration, i. e., tbe area tj, fs bounded by the lines which represent 
the four roots to the radicals In the integrand. Notice also that the coZue of UudauUe integnl 
Is ImUpeiutent of r and s and depends only on the parameter ^ the slope of the leading edges 
in tha r,i plane. 
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2E 

1 -(-mjS 


(140) 


where E has the modulus — and returning to 

the Cai-tesian coordinate system, one finds for the values 
OE ^ and <f>i, 


<pi 


4Voa 

T(l+m/3) 




(141) 




4 V^a 
\(l+m0) 




4E ■ 1 
r(l+m(3)J 


(142) 


The process of iteration could be continued, and <pu would 
be expressed as an infinite series of teims containing the 
parameter m/S. However, since the teims in the series 
expansion are all independent of x and y, or, in the character- 
istic system, of r and s, it is more efficient to write </> in the 
form 


and determine the magnitude of A by substituting this 
expression into equation (135). There results the equality 

._ 4 7og A / 4tE A 

fl-Mo + V 

from which it can be shown that 

Finally, therefore, the velocity potential on the upper 
surface of a triangular wing with subsonic leading edges 
can be wTitten 


<Pu(x,y)=^-^Jm*x^-y* (145) 

and the familiar expression for the loading coefficient (see, 
e. g., reference 10) follows immediately 

A® 4:am^z ,, , 

The purpose of examining this particular problem was 
not, of course, to obtain the solution presented as equation 
(145) or (146), since that solution is by now quite well 
known. Rather, the purpose was to show how the super- 
sonic doublet equation could be transformed to a singular- 
integral equation of the second kind and how this equation 
could, in turn, be solved by applying an iteration process. 
Such a method has far more general applications than are 
given here and is by no means limited to problems in which 
the flow is conical or quasi-conical. 

PART IV— APPLICATIONS 
DIRECT PROBLEMS 

The following tluree examples will serve to illustrate how 
the formulas derived in Part I can be used to solve the direct 
problems outlined in Part II. The applications will be 
limited to supersonic flow problems. 


RECTANGULAR WING WITH BICONVEX SECTION 

Consider a rectangular, nonlifting wing with a chordwise 
section given everywhere by the equation 

^J,x,y)=^{c-2x) (147) 

c 

where X« is the slope of the upper suirface, c is the chord, and 
t is the maximum thickness. The equation for the pressure 
on the surface of such a wing will change form in each of the 
four regions indicated in figure 16. In region 1 the pressure 
is the same as for a two-dimensional wing with the same sec.- 
tion, and the remaining regions contain the three-dimensional 
or tip effects. 


Vo 

i 



Seciion j rt 



Fiouez 16.— Rectangular nonlJfUng wing witli biconvex section, 

The equation for pressure coefficient on the wing can be 
determined from either equation (76a) or (76b). Equation 
(76a), for example, becomes for region 1 

From the result given as equation (17), it is apparent that 
the integral term is zero and the pressure coefficient in region 
1 of figure 16 is simply 

C^~,{c-2x) (148) 

One can easily show that equation (76b) yields the same 
result. 

The evaluation of pressure coefficient in regions 2,3, ami 4 
can be canied out in a similar fashion. A slightly difTeront 
approach can be used, however, that is useful in obtaining 
results in this and similar problems. Consider a wing ivith 
regions as shown in figure 17. Regions a, b, and c include 
all the area ahead of the line Xi=x. It is obvious that ro* is 
a pure real number for all Xi, yi inside region b (i. e., inside 
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Vfl 



the forecone from x, y) and a pure imaginary number for all 
ail, yi inside regions a and c. Further, it is clear from physical 
considerations that is always real on the wing plan form. 
Hence, if the symbol n is introduced to represent the com- 
bined area in regions a, b, and c, the area r in equations (76) 


can be replaced by the area Tf and the real part of the result 
will be the correct answer for the pressure coefficient.® 

By means of this concept, pressure coefficient for aU 
points on the wing can be written in the form 


(c— 2ii) (r— Xi) 


C,(x,y)=^, (c-2x)-^ R.P. dy^ 
or, alternatively, 


(149) 

(150) 


where the letters R. P. indicate that the real part of the inte- 
gral is to be taken. Evaluating, for example, the former 
equation one finds 




— ,R.P. - 

vc* L 


e—2x 


|3 


arc cos 


P(s—y) 


2(jr-s) arc cosh 

c— 2r P(8+y) I \ 1 . a: "] 

__ arc cos _^_2(y-Hs) arc cosh^^J 

(151) 

The real and imaginary parts of the arc cosine and arc cosh 
terms in this equation are given in the following table for 
all real values of the argument. 

>Tli» areas d and ecotild also be fodudedbi the definlttoii of rt since n^Iil tmBgfnaryffttbtn 
these regions also.- Hovrever, region f most not be Included since there the term nf is again 
real The definlUon of rj adopted Is nsnally the moirt craiTEtdent. 


Range of X 

— l>as 


1<Z 

arc cosh x= 

U dt 

Ji v«*-i 

tV+arc cosh. ( — a:) 

i arc cos x 

arc cosh x 

arc cos ®=J 

'I dt 

X— t are cosh (—a:) 

arc cos x 

i arc cosh x ■ 

are sin J 

^ dt 


arc sia x 

3 ,,.ccosh. 

Vl-P 




Equation (151) contains, at once, the entire solution to the 
wing shown in figure 16. Thus, in region 1 none of the 
terms in the braces has a real part. Hence, equation (148) 
follows immediately. In region 3 the solution can be 
written 

2^(y+s) arc cosh 

The solution in region 2 follows from the one given in region 
3 by symmetry; and, finally, the expre^ion for the pressure 
coefficient in region 4 is given by equation (151) wherein 
every term is real. 

DRAG KEVEBSIBUrTT THEOSEhl 

The well-known theorem that the drag of a symmetrical 

sharp-edged, nonlifting body is the same in forward and 


reversed flight at the same speed (see references 11 and 121 
can be derived in another way using the methods described 
in the above sections. 

By definition, drag coefficient is 

JJ 2X„(r, y) Op{x, y) dx dy (152) 

where S is the area of the wing. Using the real part concept 
outlined in the discussion of the preceding example, one can 
write 

JJ X„(r, y) dx dy R.P. dyi da:, 

(153) 

The equation for the drag coefficient in reversed flight can be 
obtained by: 

1. Replacing the area n by th where Tf4-Tff=;S' 
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2. Eotating tlie axial system in tlio.a^ plane tlirough 180’ 
;3. Eeversing the signs of Xu(a:, y) and Xu(a:i, y^) 

There results 


Va 



^ ^K{x,y)dxdy'R.7. ^~dxi 


(z,—x)K(xi, yO 




n54) 


Subtracting equation (154) fi-om (153) gives 

J'i!, J * J f 

,yi)K(x,y) 




(155) 


Since the symbols Xi, y^, x, y are dummj’ variables of integra- 
tion, the last teiTn in equation (155) can be written 


4 

vS 


E.P. |* J* dy 


x\^{x, y)\u{xi, yi) 


But reversing tlie operators fdy^Sdxi and J'dyJ'dx (always 
preserving the same order within the operation) and sub- 
tracting gives 

since the residual is zero. Hence, the second teim in equation 
(155) is the same expression as the first except for the sign and 

Cd — C^i>,=0 
or 

(156) 


as was to be sliown. 


LIFT ON WINGS WITH SUPERSONIC EDGES 

The lift on any wing can be Avritten. . 

(167) 

Moreover, 


where T. E. and L. E. denote the trailing and leading edges, 
respectively, and vt.e. is the value of the velocity poten- 
tial on tlie upper surface of the w'ing at the trailing edge. 

Consider now a Aving with all edges supersonic and a 
straight trailing edge not uece^arily at right angles to the 
free-stream direction (see fig. 18). Let the wing be a plate 
having arbitrary twist and camber. Then, for a point on 
the wing, the velocity potential can be widtten in the x, y 
coordinate system on the basis of equation (121) as 


X 


L.E. 2 


_4<PT.-i 

" Ho 




,^R.P. f (158) 

IT Jri J V(a:— ari)'— i3%— yO* 



Figure 18.— Wing with supersonic edges. 

where T{, as in the previous examples, is the area on the 
wing ahead of the line Xi=x. If the equation of the trailing 
edge is 

Xi=a+yi tan A 

where fl- is some constant, the value of the potential at the 
trailing edge can be written 


(PT.E.“^"" E. 


o ' 


ixuyddXjdyi 


(a+y tan A— a:,)‘-'/S%— yd’ 


(159) 


where the area of integration is the whole wing plan form 
since the trailing edge of the wing is supersonic and the 
aftercone from the point at which e. is being evaluated 
caimot intersect the whig. The total lift L on the wing ran, 
therefore, be Avritten in the form 


-=^E.P, 
2 »’Ho 


o ’ 


w„(^xj,yi)dxidyi 
V+y tan a'— ard*— |8*(y— y,)* 


( 100 ) 

The area S docs not depend on y, so tlie y integiation can 
be made fii’st and, since the edges are all supersonic, the 
interval Si<y <«2 must ahvajs contain the roots Xi and Xj 
of the expression under the radical. Hence 


E.P 


f 'tt 

tan’. 


dy 


•1 VOS'*— tan’A)(Xi-y)(y-Xd 
dy 


r__ 

J*-! V03‘'-tan®A) (Xi— y) (y— Xd 

and since 


then 


J>‘i V(Xi— y)Ty— Xd 


-4 

2 Vl3*-tan®A 


( 101 ) 


Defining the average angle of attack a bj' the expression 




—Wu{xuy\)dxidyi 

Ho 


( 102 ) 
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one can. Avrite equation (161) in the alternative form 


Cn 


4a 

Vi8*-tan*A 


(163) 


It is interesting to notice that the lift coefficient for the 
wing just studied is the same as that for a two-dimensional 
Sat plate Syiug at an angle of attack a into a free stream, 
the speed of which is given by the component of velocity 
normal to the trailing edge of the three-dimensional wing. 
This result has been derived previously in reference 13. 

INVERSE PROBLEMS 
LOW-ASPECr-EATIO EECTANGULAH WING 

It win be noted in the summary of results for the funda- 
mental formulas applicable to supersonic Sow (equations 
(75) through (82)) that the results are presented for both 
orders of integration in the x, y coordinate system. While 
the analysis of direct problems can be carried out in all 
cases if, say, the Xi integration is always performed first, it 
may sometimes be more convenient to perform the yi inte- 
vgration-first. In. the iinalysis of. inverse -.problems, however, 
it is much more important that freedom exists in the choice 
of the first variable of integration. A good example of this 
is provided by the foUowing approximate derivation -of the 
loading on a slender (in the streamwise sense) rectangular 
flat plate. 

By considering the special case when the wing chord is 
long compared to the span, one can obtain an approximate 
solution for the slender flat plate by assuming the loading 
coefficient has the form 





Figtjse 19 .— Integratiod area for slender rectangnfar wlnf . 

Introduce the notation 

. S S a — 01 



(164) 


and these equations become,® since a=— U'u/T’i), 
for O<0<1 


where a is the semispan (see fig. 19) and/(a;/s) is an unknown 
fimction. The function / is to be determined by the con- 
dition that is constant along the centerline of the wing. 
If the solution to such a problem is to be determined by Use 
of the doublet or vortex equation, it is obviously important 
that the first integration be made with respect to yi since 
the variation of Ap/g with i/i is knoAvn. 

Since Ap/g is given, let the vortex equation be used and 
let equation (164) be placed into (77a). For y=0 (and for 
added simplicity for i3=l) the area t is shoAvn by the shaded 
area in figure 19 and the resulting equation can be vrritten 
for a;>s 






yi— (Pi/a)® 


for 0<(i<Cs 





X"' yi— (3/i/g)^ 

T- cx-x,) — aji)*— yy 


dyi 


1 =/(0) -t-- r de; (1 6f.a) 

Wjo 

for 1<C6 

i=/(0)+2r knBj{e;)de^+^ r~"' EjiBdde, (lesb) 
wjs-i irjo 

Equations (165) are integral equations of the second kind 
(more specifically, Volterra's mtegral equations of the second 
kind) and the kernels are regular and bounded everywhere 
in the interval of integration. Hence, their solution can be 
determined readily by numerical processes. This has been 
done and the result in terms of the loading coefficient on 

the center line ^from equation (164) 

is shown in figure 20. 

For the purpose of comparison, the exact linearized value 
also is shoAvn in the interval where it is knoAvn, together 
Avith another approximate solution obtained by SteAvartson 
(reference 14) using a different approach. Near the leading 
edge, where the comparison Avith the exact results can be 
made, the agreement between the exact and approximate 
solutions obtained herein wfll be poorest because in this 
region the spanwuse variation deviates most radically from 
the value assumed in the construction of the integral equation. 


' The symbols B and E Indicate rillptic fiitegrsls^ See the table of symbols. 
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CONICAL AND QUASI-CONICAL FLOW 

Equation (108) and its solution, equation (111), occur 
repeatedly in the study of aerodynamics. In fact, equation 
(108) is often referred to as the airfoil equation since it 
plays a dominant role in the development of linearized, 
two-dimensional, subsonic wing theory. It appears also in 
the study of slender wings (reference 15) or wings flying at 
near sonic speeds (reference 16) since the boundary condi- 
tions lead again to the required inversion of the same type 
of integral relationship. In the present section, problems 
arising in supersonic conical or quasi-conical flow fields will 
be reduced also to this basic equation. 

Several methods exist whereby the solution to conical flow 
problems can be determined. The one to be studied here is 
based on the construction of conical elements extending 
radially from the apex of the field and inclined at an angle 
arc tan m to the x axis (see fig. 21(a)). In order to obtain 
such an element it is sufficient to subtract two plan forms 
of prescribed loading or thickness, each plan form having one 
side directed along the a: axis while the other sides are inclined 
at angles that differ only infinitesimally. 

Consider first the construction of a quasi-conical,''’ radial, 
lifting element that carries a load given by the expression 

(166) 

where C is a constant. The upwash field of a triangular 
plan form such as the one shown in figure 21 (b) can be found 
by integrating elementary horseshoe vortices over the appro- 
priate area t. Thus, using equation (77b), performing an 
integration, and making the substitutions 


0=mfi, v=^yfx, Vi=Pyi/x (167) 


Vo 


V 




(a) Lilting element. 

Cb) Triangular wing. 

Ficueb 21.— Construction ol conical cldment. 


Equations (168) can be written in the functional notation 
Wu=M ij) 

It follows that if the analj'sis were repeated for a wing with 
a slightly larger apex angle, there would have resulted 

Wu=f(6-\-A6, ri) 


one finds the result 


\x) Vo 4ir X 

. l-T 

4tX 


iv 


Vi‘dvi j( _^Y 

rj-mySlV e) 


J7 

iv 


i'‘dvi If. tjiV. 


(tj — 7?l)”,O<0, 1J<0 
(168a) 

(v— ’7i)S0<e, e<.v 

(168b) 


Subtracting these two expressions for gives the increment 
in vertical induced velocity due to a quasi-conical element 

dw»=lim A8=M do (169) 

tfl Ac/ Ou 

Carrying out the operation indicated by equation (109), 
makin g the substitution 


10 ^ conical field la defined aa one fn whlcli the Indaced velocities are ccaiatant along rays 
through a point. In the aubseciaentaiialyala this oorresponds to the case when Qoaal* 
oonlcal fields are those in whidi ic la greater than zero. 
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and distributing these elements between 6 i and flo with 
weight G( 6 ), one can finally show that the equation 


//sy uv 

-)3(«+i) r 

Uy' Vo 

4ir Jg 


jy>C(d)H(8,7,)de (170) 


applies where 




The function H( 8 , ri) has a simple pole at 6 =ri and the 
integral expression for is therefore evaluated as a Cauchy 
principal part. 

The boundary' condition to be satisfied by equation (170) 
is that T”o is a given polynomial of degree k in the 

variable 17. Hence, equation (170) is a singular integral 
equation with a pole of the same order as that in the airfoil 
equation. In its present form the equation appears some- 
what formidable, but it can be simplified considerably by a 
simple operation. Since ^'^wJx'^Vo is a polynomial of degree 
K, it follows that the (k-[- 1)®^ derivative of the right-hand 
member of equation (170) must vanish. Thus, iising the 
concept of the generalized principal part, one has 


8 ^C(^ 

IJ Jh 


a 71) 


which, by definition, can be put in the form 


J«x 


The fimction 8 ‘C( 8 ) is therefore to be found through the 
inversion of the integral equation 


X 


*0 8^0 (8) d8 


8—7} 



(172) 


But equation (172) is precisely the airfoil equation and its 
solution is given by equation (111). Hence, 


8^0 (fl) = — ■ . 

77 ^- 7 / ( 60 - 0){d-8d 

so that, finally, 



^ UfTi S'(®o— ’i) (.V — fli) 
0—17 



(173) 


Ap 1 ' X Y 6,9* 

2 raV0o— fi)(ff-0i) 


(174) 


where the coefficients bi are functions of the constants 6 ^ and 
81 but not of 6 . These coefficients must be determined from 
known conditions about the surface geometry. 

Consider the imyawed e., 81 = — 80 ) - triangular wing 
shown in figure 22. If the loading is to be determined on a 
flat plate with such a plan form ^^wJx^Vo becomes —a, 
and equation (174) reduces to 



u»ii/ Vg =-ax 

nju/Vo— ay®/S* 



Ca) Wing plan Conn. 

Cb) LoaddlstxlbatlDnBtz»l. 

FloiniE 22.— Load, dlstrlbatton on triangular wtngarfth spedQed twbt. 


Ap bp+biB 
2 


(175) 


By symmetry, the coefficient 61 in equation (175) must be 
zero and bo can be evaluated by placing equation (175) into 
(170) and integrating. There results the solution for the 
pressure coeffident shown in figure 22 where a=a and already 
presented as equation (146). 

One can go on to show that if, as for a wing pitching about 
the apex at a rate Q, (see reference 18) , 


w, 


then 


_ n-r ~^Q 

x‘Vo~ Vo ’ 


(176) 


Ap 4a:Q 28q^-6P 

2 , l-20o" ~| ^ 

_l- 0 o" 1 - 00 * _ 


where the co mplete e lliptic int^aJs K and Shave a modulus 
equal to i=Vl — 9 q*- Burther, if 


then 


‘!^=hv^ K — 2 

Vo ^'x‘Vo 


(178) 




WK-0.+eo’^E]+[-{S-eo^K+{4.-28o^E]8^ _ | 
[-50a*s:'+80o"(l +0o*)s:s-f (40o*-19eo‘'+4)S*] 

( 179 ) 
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w'here k=-jjl — 6a^ (=Vl— is the modulus of the com- 
plete elliptic integrals. These solutions are all shown in 
figure 22, where in one case a=QjVa and in the other o= 

If the function jS^Wu/aj'Fo is discontinuous but a polynomial 
in each interval of continuity, the solution given as equation 
(173) still applies. For example, consider the. case when 
(£=0 and Wv/Vo is a constant that changes sign in cross- 
ing the a: axis. (See fig. 23.) For simplicity let 6i=—6o, then 
by equation (173) 


«=0 (180) 

F 0 |i7| 




Figubk 23.— Loading on difTerentlallr deflected triangular plate. 


Again, as in the development of equation (175), 

Ap^ bg+bid 

Q. 


Now, however, the solution must exhibit odd symmetrj" and 
the constant is zero. The constant bi can be evaluated by 
substituting equation (181) into equation (170). There 
results, finally, 

(182) 

This solution is shovm graphically in figure 23, 

The methods described above can also be applied to prob- 
lems involving wings •with thickness and without lift. In these 
cases one constructs a radial element emanating from the 
origin and possessing a quasi-conical thickness distribution 




(183) 


The derivation of the iuduced pressure field associated with 
the element follows closely the analysis in the lifting case. 
First a triangular plan form is considered W'ith one side hav- 
ing the slope m and the other parallel to the frce-stream 
direction. The thickness is assumed to have the form given 
by equation (183) so that the pressure coefficient can be 
obtained from the equation 

Make the same notational changes as in the analysis of the 
lifting case; construct the element by taking the partial de- 
rivative with respect to fl; and, finally, distribute WTughted 
elementfi_DYer the w'ing plan form by making C a function of 
6 and integrating with respect to d between the limits 9, 
and 00- _ There results 


where 


iY n _ -2Q^+l ) 


e^Ciff)n^(9,v)dB ( 185 ) 




The boundary conditions require (j3/a:)‘C7, to be a poly 
nomial of degree xin ij. If the (x-1-1)** derivative of equa' 
tion (185) is set equal to zero, the relation 


_C“o 0*(7(0)de /&Y+1 1 r>«e‘Cme 

Jh (x+l)I A ' ff-v 


(18G) 


results. The function C(0) satisfies the same integral equa- 
tion that arose in the lifting case, and the solution can there- 
fore be ■w'ritten immediately in the form 



bt6‘ 


(187) 


Again the coefficients bt must be determined from known 
conditions about the surface geometry. 

Consider first the case when the pressure is constant over 
an unyawed triangidar plan foim as in figure 24. Thus 




and (b, being zero by sjunmetry) 


dzu _ bp 

d^ 


(188) 

(189) 


Evaluating the constant Jo by substituting equation (189) 
into (185), it can be sho\ni that the surface ordinate is 




(l-<?o*)C, 




2m\K-E) 


ylrn}x*—y'* 


(190) 


where the modulus of K and E is — This rcsult;- 
which is the equation of an elliptical section, is shown in 
figure 24^ 
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fcj 


(a) Triiingnlnr \f lag Tlthsiibsople edges. 

(b] Section. AAforCrCOQstant. 

Co] Section AA. for Ct TuryinsIIneaEly nrltlix. 

PiocBE 24.— Symmetrical triangniar wing with specified pressure distrlbutfons. 


Knallj, consider an unvawed triangular vting for which 
the pressure varies linearly in the x direction. For such a 
ease 


I 


Cp=(Op^x, (c=l 
and (6i again being zero by symmetry') - 

dzu X 6o+&2fl‘ 

from which it immediately follows by integration 


(191) 

(192) 




i ibo+2^-m-b^ arc cosh ■ 


Placing equation (192) into (185) and integrating, one can 
eventually show 


{ h[2K-E(5-e,^]+b,[ieo^+e,^K-2eo^E\ } 

(194) 

It is immediately apparent that the wing shape required to 
support a linear pressure gradient in the x direction is not 
unique, that there are, in fact, an infinite number, of shapes 
that will induce the same pressure distribution. (The con- 
verse, however, is not true. That is, a given shape has only 
one possible distribution of pressure.) 

Squire (reference 17) considered the thickness distribution 
that is obtained by neglecting the arc hyperbolic fimction in 
equation (193). BQs result corresponds to the case when &o 
is — 2j3*m^5a and can be written specifically 

y* (lS5a) 


^?Po=^5(r^K3-0o^£:-p;(4-20o">] (19Sb) 

where 6o/j3* can be related to the thickness chord ratio of the 
wing and the apex angle. Figure 24 shows how equation 
(194) can be used to obtain several triangular wing shapes 
all having the same linear pressiue distribution. (It is 
interesting to note that no combination of 6a and 6« exists 
that will give a real wing shape with zero pressure coefficient 
since the resulting negative ordinates would require the 
surface to cross itself.) 


AiiEs Aeroxautic.a.l Laboratory, 

National Advisory Coixmittee for Aerox.a.utics, 
iloFFETT Field, Calif., Oct. 16, 1950. 
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